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Seismographs operated successfully at the Apollo 12, 14, 15, and 
16 landing sites ror a duration between 5.5 and 8 years. Although the 
study of the seismic data revealed more about the lunar interior than 
other data produced by the Aprilo Scientific Experiments, the sparse 
lunar seismic network, the weak seismic sources, the strong scattering 
of the seismic waves, and occasional noise bursts impose severe limita- 
tions on enlarging our knowledge about the velocity distributions in the 
moon. ' In order to learn more about the physical state of the moon, much 
effort has been spent to improve the techniques for identifying direct 
shear and secondary wave arrivals on the lunar seismograms. In this 
study, inverse filters were designed to correct for the instrumental 
response and for near-surface effects and, thus, to facilitate the com- 
parison of the ground motion at the various sites. 

The constants of the long-period lunar seismographs were deter- 
mined from a least-squares inversion of the observed calibration pulses. 
The noise, especially the digitizing noise, was reduced by using two 
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positive and two negative calibration pulses. Knowledge of the seismo- 
graph parameters permits calculation of the transfer functions for 
each seismograph and accurate correction of the seismograms for the 
instrumental response. 

The Instrumental deconvolution was carried out In the frequency 
domain after removal of the predictable digitizing errors from the 
seismograms. The amplitude spectra of the Instrument-deconvolved seismic 
records revealed that other important filtering effects characteristic 
of each recording site remained in the lunar seismograms. In order to 
remove these effects, whitening filters were designed by averaging the 
spectral amplitudes of the largest long-period seismic signals component 
by component. These average amplitude spectra where shown to be 
primarily affected by the coupling of the seismometer to the ground and 
by the near-surface structure. The relative amplifications between the 
horizontal components were also estimated from these average amplitude 
spectra. 

Average horizontal -to-verti cal spectral amplitude ratios 
computed for a number of major seismic impact signals were compared with 
spectral ratios calculated for fundamental -mode Rayleigh waves in media 
consisting of homogeneous, isotropic, horizontal layers. The shear 
velocities of the best fitting models at the different sites resemble 
each other and differ from the average for all sites by not more than 
20%, except for the bottom layer at station 14. The shear velocities 
at the various sites increase from 40 m/sec at the surface to about 
400 m/sec at depths between 95 and 160 meters. Within this depth range, 
the velocity-depth functions are well represented by two piecewise linear 
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segments; although the presence of first order discontinuities cannot be 
ruled out. 

The seismograms corrected for both Instrumental response and 
near-surface effects were found to be comparable at the various sites 
and were used to obtain readings for the direct shear and secondary wave 
arrivals. To explore the constraints that direct arrival times could 
Impose on the velocity distributions In the lunar Interior, the largest 
natural impacts and shallow moonquakes were located by two velocity 
models. The first model consisted of a thicker crust (55 km) and higher 
congressional and shear velocities (8.1 and 4.6 km/sec, respectively) In 
the upper mantle, and the second model had a thinner crust (45 km) and 
lower upper-mantle velocities (7.7 and 4.4. km/sec). The arrival times 
observed for Impacts and shallow moonquakes could not be used to distin- 
guish between these two models because the travel time residuals for 
both models are similar and are well within the uncertainty of the 
observed travel times. Although the identification of the secondary 
arrivals suggests that a velocity discontinuity exists at a depth of 
about 25 km and that the crust is thinner than 55 km, no conclusive 
evidence was found for a discontinuity in the mantle or for laterally 
varying crustal thickness. 
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INTRODUCTION 


1. Research objectives. 

Since most of the information related to the Interior of the Earth 
was obtained either directly or Indirectly from the study of the travel 
times, amplitudes, and apparent velocities of seismic body waves and dis- 
persion of seismic surface waves, it was anticipated that the Apollo 
passive seismic experiments would provide the seismic velocity distribu- 
tion in the lunar interior. These velocities as well as other constraints 
were expected to contribute to the development of models of the composition 
and temperature distribution in the moon. 

Although the first decade of lunar seismology greatly expanded our 
knowledge concerning the physical state of the interior of the moon, the 
primary goal of studying the lunar seismograms remained the same. De- 
tailed studies of travel times and amplitudes lead to more and more re- 
fined velocity distributions in the moon. Owing to the small number of 
seismic stations (4) in the Apollo network, these studies required appli- 
cation of some new techniques in order to extract more information from 
the lunar seismic signals. Improvement in techniques for identification 
of shear and secondary arrivals, as well as maximum use of amplitude 
information were obvious needs. -This dissertation addressess the problem 
of designing inverse filters to correct for the effects of instrumental 
response, coupling of the seismometer to the ground, and near-surface 
structure, thus making the comparison among the various components of 
the ground motion at the four lunar seismometer sites possible. This, 
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as well as conventional filtering techniques , applied In order to decrease 
the noise contamination of the seismic signals, should result In better 
Identification of seismic phases. 

The steps to achieve these goals were as follows. 

1. Select the most energetic events from the complete data set 

2. Obtain Instrumental parameters and transfer functions 

3. Remove predictable digitizing errors from the seismograms 

4. Correct seismic records for Instrumental response 

5. Design whitening filters to remove near-surface effects 

6. Apply these whitening filters 

7. Determine the compressional and shear wave arrival times from the 
original and deconvolved seismograms 

8. Locate the seismic events with these seismic wave arrivals by 
assuming the compressional and shear velocity distributions in the moon 

9. Apply conventional filters to enhance seismic arrivals 

10. Make seismic sections and obtain consistent secondary arrivals 

11. Refine velocity distributions and determine discontinuities 

12. Repeat steps 8 through 11 with the Improved velocity models 

Although the primary objective is to obtain secondary arrivals 
and Interpret them with respect to the Internal structure of the moon, 
the present analysis may yield other results: 

a. ) Quantitative description of the differences among the seismometers 

b. ) Estimation of noise; determination of the frequency range in which one 

may expect to find seismic wave arrivals 

c. ) Empirical amplitude ratios of the horizontal components of the 

ground motion at a given station 


3 


d. ) Structural parameters of the near-surface zone 

e. ) More precise locations of the seismic events 

f. ) Station corrections and quantitative description of lateral hetero- 

geneities 

The dissertation is divided into four parts. Part 1 describes 
how the instrumental parameters and transfer function were obtained. 

Part 2 deals with inverse filtering, estimation of noise, and designing 
the whitening filters. The average horizontal -to-verti cal spectral 
amplitude ratios are interpreted In terms of structural properties of the 
near-surface zone in Part 3. The Improved arrival times for the various 
seismic phases and their implication for the deeper lunar interior are 
discussed in Part 4. The following section briefly reviews existing 
studies in lunar seismology. 

2. Literature review. 

Seismic stations were deployed by the astronauts at the Apollo 11, 
12, 14, 15, and 16 landing sites. While station 11 operated for several 
weeks and station 12 for almost eight years, stations 12, 14, 15, and 16 
were operating simultaneously for about five and a half years. This four 
station network spans the near face of the moon in an approximately equi- 
lateral triangle with a separation of about 1100 km between the corners 
with stations 12 and 14 occupying one of the corners. Each station has 
one short-period and three long-period seismometers. Table 1 gives the 
location, orientation, and operational history of these seismometers 
(e.g., Lauderdale and Eichelman, 1974). Further details on the instru- 
mentation are given in Part 1. 
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Most of the features of t!*,« lunar seismograms were totally 
unexpected by seismologists, who found these slesmograms quite differ- 
ent from those of earthquake seismograms. The long duration, contin- 
uous reverberation, emergent body wave beginnings, the almost complete 
lack of correlation between various displacement components, and the lack 
of coherently dispersed seismic surface waves were explained by relative- 
ly small source magnitudes. Intensive scattering, rapidly Increasing 
velocities with depth near the surface, and very low seismic wave attenu- 
ation (e.g.. Gold and Soter, 1970; Latham et ah , 1969 and 1970). The 
latter three features are related to the lack of volatiles In the 
heterogeneous, porous material which makes up the outer several kilo- 
meters of the moon. An uppermost limit of 20 km is given for the 
thickness of this scattering zone by Toksoz et al_. (1972). 

Rayleigh's perturbation theory (Rayleigh, 1945, p. 149) as used 
by Miles (1960), diffusion theory as used by Wesley (1965), by Latham 
et al_. (1970), and by Warren (1972), as well as experiments by Dainty 
et aU (1974b) showed that the raypaths of seismic wave arrivals remain 
unchanged from those predicted by ray theory, although the amplitudes of 
these arrivals are diminished greatly by intensive scattering. Thus, it 
was concluded (Toksoz et aU, 1974) that conventional travel time and 
amplitude studies were useful starting points in determining the velocity 
distribution of the Moon. This work was aided greatly by the availabili- 
ty of artificial Impacts of known location and origin time as well as of 
natural seismic sources. 

Classification of natural lunar seismic sources is based upon 
characteristic envelopes of the whole seismic signal (see Figure 1), upon 
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the relative strength of the shear wave, and upon the waveform matching 
of particular signals. Nearly Identical seismograms which are repeatedly 
observed at a given site must be produced by the same source mechanism 
repeating at the same site. Two types of such matching sources were 
found: thermal moonquakes and deep moonquakes. Almost all of the natural 
events can be classified Into these two types plus meteoroid Impact and 
shallow moonquake categories. 

The thermal moonquakes can be Identified by matching the short- 
period waveforms and by 29.5-day periodicity of occurence (Duennebler and 
Sutton, 1974; Cooper and Kovach, 1975; Duennebler, 1976). Each of these 
signals Is recorded at one station only. The source of these signals is 
believed to be dislocation caused by thermal stress on, or very close to, 
the surface at distances not greater than several kilometers from the re- 
cording site. Most of the repeating waveforms did not change In four to 
five years, and the source locations are associated to some degree with 
large rocks and craters (Duennebler, 1976). Thus, It appears that the 
source of thermal moonquakes Is propagation of small cracks In exposed 
rocks and, possibly, slumping of material along steep slopes (ibid.) 

The deep moonquakes are the most numerous type of seismic sources 
In the moon (Latham et aU , 1972 and 1973; Lammleln et. ah , 1974; 
tammlein, 1977; Nakamura, 1978). More than 100 different deep moonquake 
categories were confirmed by their nearly identical long-period waveforms. 
These signals are less affected by Intensive scattering than the meteroid 
impact signals and exhibit the best shear wave arrivals noticeable on the 
lunar seismograms. The deep moonquakes occur at depths between 700 and 
1100 km, below which partial melting possibly begins (Nakamura et al.,1373). 
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The times of occurence of these deep moonquakes correlate very well with 
the variation of the tidal forces In the moon, Nakamura (1978) demc n - 
strated that the energy released by deep moonquakes Is controlled by the 
tidal stress field and does not require tectonically accumulated strain 
energy. The matching waveforms within each source group enabled investi- 
gators to Improve the signal -to-nolse ratio of these signals by simple 
stacking (Nakamura et ah, 1976; Lammleln, 1977; Goins, 1978; Goins et 
ah, 1978). 

The most energetic moonquakes are those which originate at shallow 
depths. The magnitude cf the largest ones Is estimated to be about 4 on 
the Gutenberg-RIchter scale (Nakamura et ah , 1974). Variation of the 
amplitude of the shallow moonquake signals with distance indicates that 
shallow moonquakes occur In the upper mantle (Nakamura et ah. 1979). Only 
four to five shallow moonquakes were found yearly. Most likely they are 
true tectonic moonquakes releasing strain energy which accumulates as the 
moon cools. 

The similarity of the envelope characteristics of the natural 
Impacts to those of artificial Impacts distinguishes Impacts from 
moonquakes. The shear wave arrivals are usually burled In the scattered 
wave train. The energies and distributions of meteoroid Impacts both In 
space and In time were studied by Latham e£ ah 0 973), Ouennebler and 
Sutton (1974b), Duennebler et ah (1975), and Dorman et.ah (1978). 

Petrologic and mlneraloglc studies, laboratory measurements of 
the returned samples, and active seismic experiments resulted In the 
determination of the velocity distributions and the composition of the 
uppermost few kilometers of the moon. The low velocities and their large 
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gradients with depth (e.g., Cooper et al.. , 1974) were attributed to the 
presence of porous* breed ated material and to the absence of volatiles 
In this zone (Gangl, 1972; Talwanl et 1973; Todd et jH.., 1973; 

Tlttman et ah* 1974). The congressional velocity Is about 100 m/sec 
at the surface and reaches 4 to 5 km/sec at depths of 1.5 to 2 km. The 
near-surface velocities and densities are described In greater detail In 
Part 3. 

The velocity distribution In the deeper crust was determined 
from the study of artificial Impact signals. Since the shear wave 
arrivals are obscured by Intensive scattering, Toksoz et aU, (1972 and 
1974) assumed the ratio of the compresslonal and shear velocities to be 
1.7. At depths greater than a few kilometers the compresslonal velocity 
Increases slowly, attaining values of about 6.1 km/sec at 20 km (Ibid). 

The lower crust, defined as the layer below this depth, has a rather 
homogeneous compresslonal velocity of about 6.8 km/sec. 

The transition between the upper and lower crust is by no means 
considered to have been uniquely defined. Possible refracted wave arrivals 
(Toksoz et a_l., 1972 and 1974) and peg-leg multiples (Goins, 1978; Goins 
et aU, 1978) indicate that there Is a first-order discontinuity at a 
depth of 20 to 25 km. This first-order discontinuity could be the conse- 
quence of compositional changes from a basaltic upper crust to gabbroic 
anorthosite In the lower crust (Mlzutani and Newbigglng, 1973; Wang et aU , 

1973). A second-order velocity discontinuity would not require composi- 
tional changes but Implies the closing of Impact- Induced microcracks In the 
basaltic crust (Todd et ah, 1973). The crust was possibly formed by 
large scale Igneous differentiation early in the history of the moon. The 
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highland soils were derived from a suite of highly feldspathlc rocks Reid, 
1974). This feldspathlc material Is believed to represent the Initial 
crust which was cratered by Impacts and was flooded by basalts In maria. 

This sequence of events Implies the existence of a first-order discontinu- 
ity above the feldspathlc layer only In mare areas (Dainty et al_. 1974). 
Goins (1978) Identified peg-leg multiples from a depth of about 20 km also 
at highland areas. If this Interpretation Is correct, the compositional 
change should not be restricted to mare areas and should be the result of 
the original crustal differentiation, and the 20-km discontinuity airo 
should have remained undisturbed by Impacts In later times. 

The thickness of the crust Is given at 55 to 60 km under stations 
12 and 14 in most studies published since the beginning of lunar seismic 
Investigation. New arrival times by Koyama and Nakamura (1979) for the 
congressional wave read from two artificial Impact records at station 12 
indicate a crustal thickness of about 45 km. Peg-leg multiples from the 
base of the crust (Dainty et aL, 1977; Goins, 1978) Imply a crustal thick- 
ness of either 60 or 90 km at station 15 and of about 75 km at station 16. 

The velocities In the lunar mantle have been determined by a pro- 
cess of successive approximation starting with a homogeneous sphere beneath 
the crust. The average velocities between depths of 60 and 920 km were de- 
termined hy Toksoz «t £1.(1974) as 8.0 and 4.2 km/sec for the congressional 
and shear waves, respectively. Lammlein et_ al_. (1974) and lammlein (1977) 
used a model of linearly decreasing velocities In the mantle 8.1 and 4.7 km/ 
sec at a depth of 60 km and of 8.0 and 4.3 km/sec at 800 km for the com- 
presslonal and shear waves, respectively. Nakamura et jH. (1974b) showed 
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that the differences between the arrival times of the shear and congres- 
sional waves require a drop In the shear wave velocity at a depth of 
about 300 km. They later refined the velocity distributions In the mantle 
by using more events and by analyzing the decay of the shear wave anvil tude 
as a function of distance (Nakamura et al.., 1976). They showed that a con- 
tinuous decrease of the velocities with depth satisfies both the travel 
time data and the shear wave amplitude data. A discontinuity at a depth 
of about 300 km was also Indicated by reflected compresslonal wave 
arrivals (Voss et ah, 1976) and by converted shear wave arrivals (Nakamura 
et ah, 1974b; Latham et ah , 1977). Indications for lateral variations 
of the seismic velocities In the upper mantle were found by Nakamura et. 
ah (1977). 

Goins (1978) and Goins et ah (1978) Inverted the observed arri- 
val times for a two-layer mantle. Compresslonal and shear velocities of 
7.7 and 4.45 km/sec were determined In the upper mantle and 7.54 and 4.25 
km/sec In the lower mantle. Although the bounds given on these velocities 
permit the velocities to be constant In the mantle, their final model 
consists of continuously decreasing velocities In the upper mantle and of 
a first-order discontinuity at a depth between 400 and 480 km. This dis- 
continuity was determined by reflected arrivals on sections of polarisa- 
tion-filtered seismic records (Dainty et ah, 1976; Goins et ah, 19/8; 

Goins, 1978). The latest results by Koyama and Nakamura (1979) also Indi- 
cate small negative velocity gradients In the entire mantle. The mantle 
velocities at a depth of 45 km were found to be 7.85 and 4.31 km/sec for 
the compresslonal and shear waves, respectively. 

Although the bounds on the velocities of the deeper lunar Interior 
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are rather large (Nakamura et ah, 1976; Goins, 1978), the various upper 
mantle velocities are consistent with a composition of olivine-pyroxene 
(Ibid.). Small negative velocity gradients can be attributed to the In- 
crease of temperature with depth and to cotnposltlonal change caused by 
Increasing pyroxene or Iron content (Nakamura et.aU, 1976). The In- 
creased attenuation, especially of shear waves, and the deep moonquake 
activity In the lower mantle support the hypothesis that the temperature 
Increase Itself can account for the velocity distributions. 

The complicated waveforms of the lunar seismic signals and low 
coherency among the different components of the ground motion make con- 
vincing Identification of reflected or converted arrivals quite difficult. 
These secondary arrivals are needed to clarify the discontinuities and to 
narrow the bounds on the velocities of the lunar Interior. The major tool 
for searching for secondary arrivals has been the particle motion filter. 
Although the application of these filters assumes that the horizontal 
seismographs are matched, Jarosch (1977) concluded that the method seems 
to work In spite of the mismatch among the Instrumental responses and of 
the small angle of Incidence for most body wave arrivals at the free surface. 
Goins (1978) scaled the signal amplitudes of the three components to the 
same level before applying the particle motion filters. Although this 
normalization may also lead to incorrectly rotated transverse aod radial 
components, Goins (1978) also concluded that It Is a reasonable approach. 


PART 1 . DETERMINATION OF THE TRANSFER FUNCTION OF THE LONG-PERIOD 
LUNAR SEISMOGRAPHS. 

1.1. Introduction. 

A seismograph measures the relative motion between the frame of 
the seismometer and the mass of the pendulum. The relation between this 
motion and the ground displacement, velocity, or acceleration can be 
described by an integral equation, a differential equation, or the 
equivalent Laplace transform equation. These equations are functions of 
frequency and of certain controlled parameters of the seismograph. 

It is often useful to obtain the actual ground motion from the 
seismic signal. This process is called instrumental deconvolution. In 
order to calculate the ground motion with reasonable accuracy, one has 
to answer two basically different groups of questions. 

1. What are the characteristic equations of the system? What 
is the accuracy of the solution to these equations under the specified 
conditions? Are the assumptions used in deriving the solution reasonable 
approximations to the actual operating conditions? What is the noise 
introduced by the instrumentation? Do significant changes occur in the 
instrumental constants during operation? 

2. What is the reliable frequency band for a meaningful 
deconvolution, that is, what is the accuracy of the deconvolution 
process? How do errors in the instrumental response and in the digi- 
tizing process propagate? What is the signal-to-noise ratio of the 
seismic records? 

Answers to the first group of questions are independent of the 
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seismic records end come from the design, testing, and calibration of 
the seismographs. Part 1 answers these questions for the long-period 
lunar seismometers. The second group of questions Is related to the 
deconvolution process Itself and Is discussed In Part 2. 

If the Laplace transform or Fourier transform equation Is 
used to solve the integral or differential equations then the 
function which describes the seismograph is called the transfer 
function whose Inverse Laplace or Fourier transform is the impulse 
response (BSth, 1973. p. 234). The transfer function or the impulse 
response can be obtained in the following ways: 

1. The empirical method Involves Fourier analysis of the 
output of the system for an input which is equivalent to a known 
ground motion and is produced by some electrical or mechanical technique. 
This output is termed the calibration pulse or signal. Shaking 

table experiments (Silverman, 1939; Kelly, 1939), sine wave simulators 
(Murphy et al_. , 1954), application of the Maxwell impedance bridge 
(Wlllmore, 1959), and tapping tests (Matumoto, 1958; MacElwane and 
Sohon, 1932) are some examples for this technique. 

2. The theoretical method consists of deriving the impulse 
response or transfer function analytically. The impulse response can 
be obtained as a solution of the equations of motion (Chakrabarty , 

1949; Chakrabarty and Choudhury, 1964; Landisman et aK , 1959). 

The transfer function may be derived by applying electrical network 
theories (Neugeberger, 1970; Kollar and Russell, 1966; Dopp, 1964; 
Hagiwara, 1958; Sutton and Latham, 1964). These methods require 
knowledge of the values of the instrumental constants. 
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The first technique does not yield the Impulse response If 
the Input Is not an Impulse In the desired domain (displacement , 
velocity* or acceleration). In this case* additional numerical 
Integration or differentiation Is needed to derive the Impulse 
response from the calibration pulse. These operations Introduce 
additional numerical noise, which would obscu^ the comparison of 
the various seismometers and would reduce the width of the frequency 
band for a reliable Instrumental deconvolution. For example. If the 
Input Is obtained from a tapping test (MacElwane and Sohon, 1932, 
p. Ill) then the impulse-displacement response of the seismograph 
Is the twice-differentiated calibration pulse. If Lagranglan inter- 
polation polynomials are used to obtain the derivatives, numerical 
noise can severely limit the passband. Increasing either the number 
of points used for the interpolation polynomials or the noise level of 
the pulse will narrow the frequency band for a predeterml ned signal-to- 
nolse ratio for the instrument-deconvolved seismic records. 

To decrease such noise and thus to broaden the reliable 
frequency band, and to learn more about the seismograph, it is desirable 
to combine the empirical and theoretical methods. First, the 
Instrumental constants are determined from the calibration pulse, and 
then these constants are used In the analytical description of the 
Impulse response of the transfer function. Espinosa et al.. (1962) 
showed that the instrumental parameters can be determined by canparing 
the calibration pulse with a set of theoretically calculated responses 
at several points only. Mitchell and Landlsman (1969) used the 
numerical least-squares technique for an electromagnetic seismograph. 
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They calculated the theoretical transfer function with a formula 
given by Haglwara (1958), and transformed It Into the time domain 
with the fast Fourier transform algorithm. Their experiments showed 
that this method was superior to the direct Fourier analysis of the 
pulse In yielding the values of the seismograph constants and In 
approximating the true transfer or response function In the presence 
of noise. The least-squares method was Improved by deriving explicit 
expressions for the theoretical Impulse response In the time domain 
(Jarosch and Curtis, 1973) and also by Including the scale factor 
between the observed and theoretical calibration pulses (Ibid.) and 
the origin time of the observed calibration pulse (Mltronovas and 
Wlelandt, 1975; Mitronovas, 1976) in the Inversion. 

In the present analysis, the least-squares technique Is used 
to determine the instrumental constants and the transfer functions of 
the direct-digitizing, feedback-controlled, long-period lunar seismo- 
graphs. These lunar seismographs are described briefly in the next 
section; Section 1.3. provides computational details. The influence 
of noise and the results of these calculations are discussed in the 
final sections. 

1.2. Description of feedback-controlled lunar seismographs. 

The Apollo lunar seismic stations have a sensor unit and a 
central station electronics module. Each sensor unit consists of one 
short-period vertical component seismometer (SPZ) with a resonant period 
of 1 second, three closely matched, orthogonal , long-period seismometers 
(LP) with resonant periods of about 15 sec, and uncaging, leveling, and 
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temperature control devices. In the short-period Instrument, a moving 
magnet mass is suspended from the frame with a LaCoste spring. By Its 
movement, this magnet Induces voltage In a coll, which is fixed to the 
frame. This voltage Is therefore proportional to the relative velocity 
between the frame and the suspended mass. A LaCoste type spring Is 
also used for the suspension of the mass In the vertical long-period 
seismometer (Z), while the horizontal sensors (X and Y) employ a 
swinging-gate system. Since the relative motion between the mass and 
the frame of the long-period seismometers Is measured by capacitor- 
type displacement transducers, the electrical output is proportional 
to the relative displacement between the mass and frame. Long-term 
stability of the Instruments Is accomplished by feedback circuits, 
which can be switched between the normal (flat) and the modified 
(peaked) modes o* operation. In the latter mode, the feedback filter 
Is bypassed and the transfer function was sharply peaked at a period 
of about 2.2 sec. 

The filtered and amplified signals of all four seismometers, as 
well as the feedback signal (also called the tidal output), and the 
output from the temperature sensor unit were digitized in 10-bit words 
and multiplexed in the central electronics unit. The 10-bit words were 
transmitted through a shift register in frames of 0.60375 seconds each. 
These frames comprise four samples from each of the LP components, 

28 (station 15 only) or 29 SPZ samples, and alternately the two 
horizontal tidal samples or the vertical tidal sample and the output 
of the temperature sensor. Although the SPZ at station 12 never worked 
and occasionally a few of the other components functioned erratically, 
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most of the sensors operated normally for approximately 5.5 to 8 years 
until the time of the turn-off of the lunar passive seismic network on 
September 30, 1977. More detailed descriptions of the Instrumentation 
are given elsewhere (e.g., Sutton and Latham, 1964; Earth Sciences, a 
Teledyne Company, 1968; Latham et,al. , 1969 and 1970; Lauderdale and 
Elchelman, 1974). 

The calibration of all twelve LP seismographs was accomplished 
periodically by commands (e.g., Latham et al., 1969 and 1970). A suddenly 
applied electrical current, I, through the damping coll at time t»t' 
produces a force, which Is equivalent to a step of ground acceleration 
(Sutton and Latham, 1964; Kollar and Russell, 1966; Jarosch and Curtis, 
1973): 


dy (t) dyj(t) dy (t) 

my a (t) * m — - — * m — =* — * ylH(t-t'), or m — - — 3 yl<5(t-t') (1A) 
a dt dr dt 


where y is the ground motion and subscripts a, v, and d denote 
acceleration, velocity, and displacement, respectively; y is the 
damping coll constant: His the Heaviside step function, s Is the 
Dirac delta or impulse function and m is the seismometer mass. The 
relative displacement measured between the mass of the long-period 
pendulum and the frame of the seismometer is the calibration pulse, 
C(t), which is proportional to the response to a unit step of 
acceleration, R Q (t). 

The response to a unit impulse of acceleration, R^(t), is 
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given by applying the differentiation and Integration theorems of 
convolution 

Rl t 

SjC(t) - R 0 (t) » H(t)*R 1 (t) = (IB) 

where * stands for convolution, x Is a dummy variable, and the system 
was assumed to be causal. From (IB) 


dR 0 (t) 

dt 


R,(t) 


OC) 


The responses to a unit impulse of velocity or unit impulse of displace- 
ment, Rg( t ) and R^(t ) » respectively, are obtained by the definition of 
the responses. 


y a (t)*R 1 (t) * y v (t)*R 2 (t) » y d (t)*R 3 (t) 


(ID) 


and by the differentiation and Integration theorems of convolution, 


dR,(t) dR,(t) d 2 R,(t) 

y a (t)*R,(t) = y,(x)dx* s y v (t)*- — = y H (t) — l — 

a 1 a dt v dt d dt^ 

(It) 

Relations among R^(t ) ' s described by equations (IB) through (IE) can be 
summarized by 


*- 0 , 1 , 2, 3 


(IF) 
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dlUt) 

or "" Sr * 11 R k+1 ^ » k-0,1 .2 (16) 

Equations (IF) and (16) can be expressed In the frequency domain as 
follows 

^jc^CU) » T k U), k*0,l ,2,3 (2A) 

ju.T k (u») - T k+1 U). k-0,1,2 (2B) 

where a) Is the angular frequency, j»/T, and (T(w) and T k (w) are the 
Fourier transforms of C(t) and R k (t), respectively, 

OD 

ITU) * f C(t)exp{-jut}dt (2C) 

OD 

T k («) » / R k (t)exp{-ju.t}dt (2D) 


If the seismograph responds linearly to the ground motion, the 
Fourier transform of the calibration pulse for the long-period lunar 
seismographs can be approximated analytically by using the equivalent 
electric current analog shown In Figure 2 (Sutton and Latham, 1964; 
Earth Sciences, a Teledyne Company, 1968) a* follows 
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jo.* 1 K 1 G(«)F d (cu) 

1+K, K^G(u)F d («)F^(u) 


<l<3F a (w)Fd(w) 


juf 1 K^G(w)Fd(w) 
l+K^GUlF^w) 


«V a (-)F d (-) 


(3A) 


(3B) 


where superscript c denotes calculated values* 


G(<u) * - y " 

£ .O * 


(3C) 


-« +26w 0 jw+u) 0 


Is the transfer function of the seismometer If the Input is considered to 
be frame acceleration and the output is the displacement between the 
seismometer mass and the frame (Sutton and Latham, 1964), 


F d<“> ■ 


(3D) 


is the transfer function of the demodulator low-pass filter, 




F f<“> * 


(3E) 


is the transfer function of the feedback low-pass filter, 


F (oj) *- 
a a) + ju) 

3 


(3F) 
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Is th« transfer function of the high-pass filter In the output 
amplifier! 

2 


Is the transfer function of the output low-pass antlallas filter, and 
K, K| , < 2 » Si Wq , and ti)| are defined In Table 2* 

Equations (3A) and (3B) are valid when the feedback circuit Is 
operated In the flat mode and the peaked mode, respectively. 

The transfer functions to an Impulse of acceleration, velocity, 
or displacement, t£(w) (k+1,2,3), can be calculated from equations 
(2) and (3) 

Tjj( w ) - (J«) k {R(»)+;!(«)} - u k A(a.)exp{j(F(u.)-kf)} (4A) 

where R(w) and I(w) are the real and Imaginary parts of the calculated 
transfer function to a step of acceleration, A(w) and F(w) are the 
calculated amplitude and phase spectra of the response to a step of 
acceleration 

TjjU) « R(u) + jl(w) - A<u)exp{jF(u)} (4B) 

A(oj) ■ \/r(u>)^ I{u») 2 , F(u>) • tan"^{I(u)/R(a»)} (4C) 





The nominal values for the -alsmograph constants used In 
equations (3) and In Figure 2 are listed In Table 2. The nominal 
responses to a step of acceleration calculated with the Inverse 
Fourier transform of equations 1(3) differ from the observed calibration 
pulses shown In Figures 3 and 4 for the peaked mode and the flat 
mode, respectively. In these figures, the positive and sign-reversed 
negative pulses were averaged. The averaging and the determination 
of the beginning of these calibration pulses are discussed In the 
following section. These pulses were taken during the relatively 
quiet lunar nights when the Instrument temperatures were 125+0. 5°F 
and the background noise was lower than that during the lunar days. 
Table 3 gives the dates when these calibration pulses were recorded 
and summarizes the visually obtained properties of these pulses. 

There are significant differences among the amplitudes of the 
pulses for the different components. Some of the differences In the 
amplitudes are attributed to the different currents applied through 
the damping colls. Although the current level was set to provide 
about one-fourth full scale output or 256 digital units (DU) In the 
flat mode (Earth Sciences, a Teledyne Company, 1968), all of the 
calibration pulses In the peaked mode, except for component 15Z, 
are much smaller (13 to 20 DU) because the effective stiffness of the 
suspension Is Increased In the peaked mode. The larger calibration 
pulse for component 15Z In both modes Is primarily due to an 
undetermined failure In setting the calibration current. The cali- 
bration pulse for component 12Z looks similar to the others in the 
peaked mode but differs from the others In the flat mode. 
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To minute the effects of temperature on the response of the 
lunar long-period seismographs and the long-term drift, the maximum 
amplitudes, durations, and time Intervals between successive peaks 
were compared for calibration pulses taken In the peaked mode between 
1971 and 1975. In the flat mode,, the times of the zero crossings 
were compared Instead of the time Intervals between successive peaks 
for calibration pulses recorded In 1975 and 1976. 

The changes In most calibration pulses for a given component 
were within the seismic background noise level except In the peaked-mode 
calibration pulses for component 15Z. For these pulses, the time Inter- 
vals between successive peaks remained unchanged and the maximum ampli- 
tudes were 14% lower at 137°F than at 125°F. Most of the differences 
can be explained by a smaller calibration current and/or output gains at 
137°F thar, at 125°F. In the following sections, those calibration pulses 
are analyzed which were taken during the lunar nights, when the tempera- 
ture was 125jp.5°F, and the signal -to-nolse ratio was highest. 

1.3. Computational details for determining the Instrumental constants 
by least-squares fit of the observed and theoretical calibration 
pulses. 

The least-squares scheme, described In the Appendix, can be 
constructed either in the time domain or in the frequency domain. 

The former scheme fits the calibration pulse with the calibrated 
response to a step of acceleration. In the 7 atter scheme, either 
the calculated amplitude or the calculated complex transfer function 
, serves as a fitting function to the Fourier transform of the observed 
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calibration pulse. The schematic flow chart of the computer program 
1$ shown In Figure 5. The first method Is preferred to the second one 
because It Involves numerical operations on the calculated transfer 
function Instead of the observed calibration pulse and because the 
observed-mlnus-calculated (0-C) curve Is measured In digital units. 
Since this scheme requires the calculation of the Inverse Fourier 
transform of the partial derivatives ♦ and thus, longer computational 
time, It was used only to Investigate the effects of noise and to 
judge the goodness of the fit. 

In order to understand the limitations of this numerical 
technique and to Interpret the results correctly, I shall 1 discuss the 
Iteratlonal process In some detail (see Figure 5 and Appendix). 

The first step In the iterational loop requires the 
calculation of the transfer function to a step of acceleration. The 
nominal values of the seismograph constants, listed In Table 2, were 
used for the first Iteration. If the sampling Interval In the time 
domain Is t(t»0. 15094 seconds for the long-period lunar seismograms) 
and the number of points to be calculated In the time domain Is N 
(M was an Integral power of 2, because the fast Fourier transform was 
used), then the Nyquist angular frequency, , and the sampling 
interval in the angular frequency domain, o, are 

■ j * 2^*3.3126 sec" 1 and a ■ ( 5 ) 

The transfer function was evaluated at angular frequencies of a ■ m c 
(m“0,l ,2,. . . ,N/*-l) according to equations (3) and (4). The analytical 
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expressions for the partial derivatives of the transfer function with 
respect to the Instrumental parameters were evaluated according to 
formulas derived from equations (3): 


3T 0 ( w) T (w) 


(6A) 


* 0 U> 

*2 

3T 0 («) 

*3 


T Q ( ci»)Gtc>>) 

0 («) 


T Q U)6U) 

D(u) 




(6B) 


(6C) 


3T U) 

V j 

jbl 

s *« 

aJUTT 

a^T 

JT 0 U) 

T („) 

m 5 

V^r 

»v«> 

2T 0 (») 

p ' re 

-2« c +j«i>2cosge-g 

3dg 

*6 

ag-u^+ju)2cosg«g 


-2 w 2 +j w 2cos| a «g 

ag-w^jo^cosg^g 


(6D) 


(6E) 


(6F) 


3T q («) T 0 («)G(w)Fj(«)F^(u) 

**7 DU) 

(66) 

3T^(w) ^(wJa^GjwjF^tu) ju 

3®g a 8 0( (it) 

(6H) 

where 0(w) ■ 1+ a 7 F^( u )F d ( u )G(u) 

(61) 
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Parameters a m (m»l,2,...,7 or 8) are given In Table 2 and F f (u) Is 
taken as unity In the peaked mode and 3T 0 (w)/3a G applies in the flat 
mode only. The real and Imaginary parts of the partial derivatives are 
shown in Figures 6 and 7. The theoretical calibration pulse and Its 
partial derivatives were calculated by the Inverse Fourier transform 
of equations (3) and (6) when the time domain was used fr- fitting 
(see Figures 8 and 9). When the amplitude spectra were fitted, the 
amplitude and phase responses to a step of acceleration were obtained 
from (3) and (4), while the partial derivatives (see Figures 10 and 11) 
were calculated with formulas derived from (4): 


3A(u) 


m 


m m 

R(o>)^ + I (u>)^ 


m=l ,2,... ,7 or 8 


(7A) 


m*l,2,...,7 or 8 (78) 

**m R(u; + l(i»i 

As seen from Figures 6 through 11, some of the partial deriva- 
tives are practically identical. The linearly dependent partial 
derivatives cannot be used in the least-squares scheme (see Appendix). 
In fact, the iteration did not converge if partial derivatives with 
respect to paramaters ag, a^, and a^ were used simultaneously in the 
peaked mode. Similarly, only one partial derivatives could be used 
from parameters a- and ag, from a^ and a g , and from a g and a 7 in the 
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flat mode. Beside the seismograph parameters, the exact origin time 
of the calibration pulse was also unknown and had to be determined 
In the least-squares scheme, which required calculation of the 
partial derivative of the calibration pulse with respect to time. 

The time derivative Is similar to the partial derivative with 
respect to a g In the peaked mode and with respect to a 4 and a g In 
the flat mode. 

The value of a^ was fixed in both modes of operation because 

its effect on the calibration pulse is negllgable. The nominal 

value of a g was also accepted as the correct one because shifting 

the beginning of the calibration pulse by one sample changes a g 

by about 25%. In the peaked mode, only the product of the damping 

constant and the free frequency of the pendulum can be determined 
2 

because is much smaller than are the other terms in the denomi- 
nator of equation (3C) at those frequencies where the amplitude 
response Is significant. I assumed 2ir/u) 0 s 15 sec in the peaked mode 
in order to determine one constant from ^ an< ^ a 3 * finally, in the 
flat mode, a^ and a 7 were adjusted to permit comparison of the 
calculated values for a 7 in the two modes of operation. Also, the 
pendulum constants are more readily affected by unusual conditions 
than are the electrical constants. 

The next step Involves the calculation of an approximate scale 
factor, c, which relates the observed points, y^ , and the theoretical 
function, y (x), such that (compare to equation (A3) in the Appendix): 
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2 

z {y< - cy(x 4 )r » minimum with respect to c only, (8) 

1»1 1 1 


where x and y are either the time and the calibration pulse or the 
frequency and the transfer function (amplitude or complex spectrum), 
and M Is the number of points to be fitted, c Is given by the 
calibration parameters (», I, and m, see equations (1A), (10), 

(3A), and (3B)): 


c 


nL 

m 


Equation (8) yields 


c 


M 

V 

1=1 

IT 

z 

1*1 


yjy(x t ) 

y(^) z 


(9) 


do) 


The sample number k closest to the origin time of the calibra- 
tion pulse was calculated by the formula 


M-k+1 

£ {C( (i+k-1 )t) 
1=1 


c 2 

cR q ( 1 t ) > =minimum with respect to k, (11) 


except when the amplitude spectra was used in the least-squares scheme. 
This criterion cannot be applied in the flat mode because the 
relative change between two consecutive amplitudes of the calibration 
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pulse Is smaller, thus the minimum expressed by (11) Is much less 
affected by the Improper beginning In the flat mode than In the peaked 
mode. 

Next, the difference between the observed and calculated 
functions, termed 0-C error or 0-C vector, and the root mean square 
(RMS) error were calculated according to equations (A2E) and (A3) of 
the Appendix. The following steps were straightforward: check for 
convergence, scale the theoretical function, and calculate the 
Inverse Fourier transform of the partial derivatives If needed. 

Finally, the adjustments were calculated according to (A3) of the 
Appendix with subroutine * LLSQ ’ of the IBM Scientific Subroutine 
Package (IBM Application Program, 1970). Before the iteration was 
repeated, the adjustments In the origin time of the calibration pulse 
was checked; if it was more than half of the sampling Interval then 
the first sample of the observed calibration pulse was changed 
accordingly. 

This Iteration process converged very fast; after 3 to 5 

-5 

interations the adjustments were of the order of 10 , and after 5 to 7 

iterations the RMS error started increasing because the limit of 
resolution was reached. 

1.4. The influence of noise. 

Special attention should be given to the effects of noise in 
the determination of the seismograph parameters because the signal -to- 
noise ratio Is rather small for most calibration pulses. Noise in the 
calibration pulse may force the least-squares scheme to adjust 
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an Instrumental parameter Incorrectly when the product of the needed 
adjustment to this parameter and the partial derivative of the spectrum 
of the calibration pulse with respect to this parameter changes with 
frequency In the same way as the noise spectrum. The following error 
sources will be discussed In some detail: 

a) Digitizing noise and Instrumental nonlinearity 

b) Instrumental noise except nonlinearity 

c) Ground noise 

d} Numerical noise 

a) The lunar seismometers used truncation -type analog-to- 
digital converters (Gold and Rader, 1969, p. 99). The digitizing errors 
can be important for the calibration pules in the peaked mode because the 
signal amplitudes were only 13 to 20 digital units (DU). Instrumental 
nonlinearity did not affect the largest calibration pulses significantly 
as demonstrated in connection with the numerical noise by the satisfact- 
ory fit obtained between the observed and calculated calibration pulses 
for station 15. The positive and the negative calibration pulses of 13 
to 20 DU are slightly different from each other in the peaked mode and 
the negative calibration pulses resulted in higher values for the damping 
constant, for the scale factor, and for 2ir/oj* than did the positive cali- 
bration pulses. From these results, I conclude that the sum of the posi- 
tive and the sign-reversed negative pulses must be used for determining 
the seismograph parameters, which also results in cancelling the d.c. 
level caused by digitization. In the flat mode, the digitizing errors 
and the nonlinear effects of the instrument are less important than 
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either the seismic ground or Instrumental noise. 

b) and c). The Instrumental noise Is described In the 
literature In great detail (e.g., Byrne, 1961; Sutton and Latham, 1964). 
fly purpose was to determine the effects of the Instrumental and seismic 
ground noises on the accuracy of the determination of the seismometer 
constants. The high frequency noise, which Is caused by the amplifier- 
transducer nols^does not affect the calculations because the instru- 
mental response Is low at high frequencies. Studies by Byrne (1961) 
and Sutton and Latham (1964) showed that the feedback and thermal 
noises increase with period and are inversely proportional to the 
square root of the mass of the pendulum. Since the mass of the lunar 
seismometers is very small (0.75 kg for the long-period instruments), 
most of the long-period noise can be attributed to these noises, 
especially to the thermal noise. This long-period noise and the 
seismic ground noise influence the determined values of a^ and, to 
a lesser extent, the scaling factor in the peaked mode (see the partial 
derivatives in Figures 8 and 10). Since the noise is about 1 DU in 
the peaked mode, the seismic ground noise cannot be separated from 
the errors of digitization. In the flat mode, the long-period ground 
noise influences all five partial derivatives used in the least-squares 
process since the corresponding constants can be determined best from 
the long-period portion of the calibration pulse spectra, where the 
amplitudes of the partial derivatives are largest. The damping 
constant is more sensitive than the other constants to noise at periods 
of 15 to 25 seconds. The rest of the parameters are influenced by 
the noise at longer periods (see Figures 9 and 11). 
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d). The numerical noises are Introduced by the finite number 
of samples used In the fast Fourier transform, by the simplification 
of the model from which the theoretical transfer functions are calcu- 
lated, and by the rounding errors. These numerical errors were 
Investigated In the time domain using calibration pulses for station 15. 
The 0-C values for the calibration pulse for the vertical component In 
the peaked mode were less than 1 DU for most of the 379 points used 
(10 points had values between 1 and 1.2, another 10 between 1.48 and 
2), If the sampling Interval In the angular frequency domain was 
ca&J Nq/4096. When the sampling Interval In the angular frequency 
domain was increased, the 0-C values Increased, but the calculated 
constants, except ag, remained essentially the same, ag increased 
by 5% or so with each doubling of the sampling interval in the 
angular frequency domain up to a=ca^/512. These observations reflect 
the sensitivity of a 5 to long-period noise (see Figures 8 and 10) 
and the long duration of the calibration pulse (Table 2). In the 
flat mode, however, the RMS errors for the horizontal components 
did not decrease for sampling interval in the angular frequency domain 
smaller than a*w^/512 and the 0-C values were comparable with the 
estimated noise in Table 3 (7.1 DU for 15X and 9.2 DU for 15Y). These 
results show that although the 0-C values are affected by the length 
of the fast Fourier transform, 1024 samples in the time domain are 
sufficient for obtaining the instrumental constants In both operating 
modes. 
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1.5. Results and discussion. 

In the previous sections, I discussed the limitations 
Inherent In calculating the constants of the long period lunar selsmo- 
graphs. They can be suimarlzed as follows: 

1. The d.c. level can be suppressed and the slgnal-to-nolse 
ratio can be Improved by averaging the positive and sign-reversed 
negative calibration pulses (I refer to this sum as a set of pulses). 

2. The transfer function and Its derivatives should have at 
least 512 complex points. 

3. The computation is faster In the frequency domain. 

4. The nominal value of 0.72 sec for 2ir/w| is acceptable as 
the best estimate because its determined value is sensitive to the 
origin time of the calibration pulse. 

5. The nominal value of 100 sec for 2 t/w. is acceptable as 

a 

the correct one in both operating modes because its derived value is 
strongly influenced by long-period noise in the peaked mode and by 
KjKg in the flat mode. 

6. In the flat mode, long-period noise affects every partial 
derivative very strongly. 

7. The scale factor, the origin time of the calibration 
pulse, the product of K-j and and the product of 8 and can be 
calculated in the peaked mode. 

8. In the flat mode, the scale factor, the origin time of the 
calibration pulse, the product of and K- which is influenced by w a 
the damping constant, and the free frequency of the pendulum which is 
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Influenced by can be adjusted. 

9. s and the scale factor are over- and underestimated using 
the negative and the positive pulses, respectively. In the peaked mode. 

Table 4 lists the average of the constants determined from the 
sum of two sets of calibration pulses (one from 1972 and another from 
1975) In the peaked mode and for one set of pulses (from 1976) In the 
flat mode. The nominal values for those constants which cannot be 
calculated are also given In Table 4. Since the determined values for 
u> o 0 and KjI <2 are different In the peaked and In the flat modes, and 
were calculated from the constants. wouTd be 1 and would be 
If the noise were negligible and the partial derivatives with 
respect to every seismograph parameter were linearly Independent. 
Relative amplifications were also calculated assuming the amplification 
to be 1.0 for component 15X, whose calibration pulse possesses a 
relatively low noise level In the flat mode (see Table 3) and whose 
determined constants are closest to the averages of the constants In 
both modes. The errors In the determination of the seismograph 
constants were calculated by taking the maximum of the differences 
between the values given in Table 4 and the values determined from one 
positive and one negative pulse only. These errors, as well as the 
relative differences In the damping coefficients, In K^, and In the 
scale factors, are given In Table 5 and correlate strongly with the 
estimated noise In the flat mode (see Table 3). The noise and the 
relative errors are larger for the horizontal components than for the 
vertical ones. 

Although the noise and the choice of the seismograph constants 
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to be adjusted In the least-squares scheme Influences the determined 
values of the seismograph constants, the transfer functions of the 
seismographs can be calculated with good accuracy. The good fit 
between the peaked-mode calculated and observed amplitude and phase 
responses to a step of acceleration is shown In Figures 12 and 13 for 
15X and 15Z. It should be noted that the relatively large 0-C values 
at long periods are caused by the finite length of the Fourier 
transform (see Section 1.3.). 

The peaked-mode amplitude responses to an Impulse of displace- 
ment, which were calculated with the constants of Table 4, with KK^-I 
DU/cm, and with equations (2)-(4), are shown In Figure 14. Tables 6 
and 7 and Figures 15 and 16 give the absolute values of the relative 
differences between these Individual amplitude responses and the 
amplitude response of components 15X and 15Z. These relative differ- 
ences are the relative errors which could be made In the deconvolution 
process by using the transfer function of 15Z, which was determined with 
the smallest errors, and of 15X, which represents the average of the 
transfer functions. As seen from Figures 15 and 16, the errors are 
most sensitive to <2 which can be determined with the best accuracy 
in the peaked mode. The maxima of the errors generally exceed 10X, which 
corresponds to the largest determined difference between the average 
and the individual damping constants (see the curve for 14Z In Figure 
15). 


1.6. Conclusions. 


Although precise seismograph constants cannot be calculated from 
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the calibration pulses* a theoretical transfer function whose Inverse 
Fourier transform fits each Individual calibration pulse best can be 
determined for each seismograph. The results allow one to deconvolve 
the records with Instrumental response having errors of less than 10X, 
out do not permit to compare the absolute amplitudes among the various 
components. Therefore, any studies using the amplitudes of the seismic 
records should be based on sound experimental and/or statistical 
Investigations of the Individual amplification factors. More precise 
determination of the transfer functions can be obtained by stacking the 
calibration pulses, by studying the long-period noises In the flat mode 
for the vertical and horizontal components separately, and by merging 
the calibration pulses of peaked and flat modes Into the same least- 
squares scheme. It Is also possible that a satisfactory explanation 
can be found for some of the Irregularities In the operation of the 
seismographs and that the temperature effects on the free period and 
damping coefficient of the seismometer can be determined by stacking 
the calibration pulses for different time intervals. 



PART 2. INSTRUMENTAL AND WHITENING DECONVOLUTION OF THE LONG-PERIOD 
LUNAR SEISMOGRAMS. 


2.1. Introduction. 

The first part of this dlssartatlon daalt with the long-period 
lunar seismographs. The theoretical transfer function of these 
Instruments was determined Individually by a least-squares process. 

This was necessary In order to correct the seismic records for 
Instrumental response In the peaked mode because the calibration 
pulses of the various components exhibit certain differences but are 
generally too small to be used directly for calculating the actual 
ground motion. The Instrumental deconvolution may increase the noise at 
those frequencies where the instrumental response is low In comparison 
to the noise. It is therefore necessary to control the noise at 
these frequencies. Different deconvolution and quality control techni- 
ques are reviewed in the following section. Then, the most suitable 
technique Is used to remove the Instrumental effects from the long- 
period lunar seismic signals recorded In the peaked mode. For signals 
recorded In the flat mode, the Influence of Instrumental response can 
be Ignored by comparison with local site effects and noise, because the 
Instrumental response Is nearly constant In the frequency band of 
Interest. Section 2.4. shows that the amplitude spectra of the ground 
motion are characteristic of each recording site. These near-surface 
effects were removed by Inverse filters derived from the average 
spectral amplitudes of the ground motion. These Inverse filters are 
termed whitening filters, because they were designed to diminish 


36 



37 

dominant peaks In the amplitude spectra. The application of these 
whitening filters permits the comparison of the various seismic 
records without the strong Instrumental and near-surface effects. 

2.2. Wavelet-deconvolution techniques. 

The seismic signals carrying Information about the transmission 
properties of the medium are modified by the source, the receiver, and 
the local structure at the recording site. The mathematical model for 
describing these modifications consists of the convolution of the 
response functions of the different linear physical processes. This 
type of model permits correction for any of these effects if Its response 
function Is given. 

The convolution equation for a linear system Is given by 
r(t) » w(t)*g(t) + n(t) (12) 

where t Is the time variable, w(t) Is the Impulse response of a system 
whose Input and output are g(t) and r(t), respectively, n(t) Is additive 
noise, and * represents the convolution operation. If w(t) Is the 
Impulse response of the seismograph, then g(t) Is the ground ration, and 
r(t) Is the seismogram. 

An Inverse filtering process can be defined as a process of 
calculating g(t) from r(t) and w(t). This process Is called wavelet 
deconvolution and can be carried out either In the time domain or In the 
frequency domain. In the time domain, the standard wavelet deconvolution 
Is performed by the Wiener or optimum least-squared filtering 
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(Robinson and Treltel, 1967). 

Tht frequency domain deconvolution can be performed after 
taking the Fourier transform of equation (12) (functions denoted by 
lower-case and capital letters are Fourier transform pairs): 

R(w) - WU)GU) + NU) (13) 

The Fourier transform of the calculated Input, ?(<n), Is obtained by 
spectral division: 

E(„) . £ial . BU l P l tt L (14) 

JTM |Bl«)r 

where superscript * denotes complex conjugate, and 7(u) is the Fourier 
transform of the wavelet used and may differ from that of the true 
wavelet, W(w): 

ST(«) • W(«) + N( u ) (15) 


or in the time domain 

w(t) • w(t) + m(t) (16) 

The error m(t) will contribute to the total error of the deconvolution. 

The relation between the Fourier transforms of the actual and 
calculated inputs Is obtained from equations (13) - (16): 

G(«) . S(u) L- Sb&MA . - p ( w )} (17A) 

1 - H(«) /lT(w) 
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where p(«) Is the relative difference between the calculated and 
actual inputs: 


p ( u ) a At «■>) “ , i _ 1“N((j)/R(tt») 

(a(u») 1-M(a»)/W(a») 


(17B) 


If M(«)*0 for all frequencies, then (17B) becomes 



The absolute values of the relative errors can be estimated from (17B) 
if |N(w) j«|R(u>) | and !M(«) |«|I(u) | : 

ip( M )l < 1 H M1 + Mali - lM(a?L + Mialt — (17D) 

' |W(u>)| |R(»)| |?( tt )| |?U)||7(«)| 

Although equation (14) is satisfactory in the noiseless case 
( |M(cu) |=0 and |N( U )|= 0 for all frequencies), |p(o>)| is large (see equa- 
tion ( 1 7D) ) at those frequencies where |W(u,)| becomes small in compari- 
son with jN(u) |/|6( w ) j (Rice, 1962). Therefore, it Is necessary to 
control the signal-to-noise ratio of the calculated input at these 
frequencies. This quality control can be achieved by suppressing either 
the nois^ NUX or the amplitudes of 1/WU). The latter process will be 
discussed in the following paragraphs. 

One can add white noise to the wavelet, w(t), before the ground 
motion is calculated. This process is called whitening and is a 
standard procedure in least-squares filtering. The equivalent process 
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In the frequency domain Is called v-Fourler filtering, which modifies 
the power spectrum In the denominator of (14) by adding a term pro- 
portional to the signal or to the noise level (Deregowskl, 1971): 


J( u ) . R(m)B*(u) V E n . ()+T)R(«)fl*(tt) 
|Sf(»)| Z «E n E s |ff( u )| Z nE s 


(18A) 


where E - ! Z du 


(188) 


E n ■ /|G(w) - G(ui) | Z dw 

■ OB 


(18C) 



It can be easily shown (ibid.) that this process increases the central 
auto-correlation value of the wavelet w(t) by the factor (1+y), just as 
the whitening does, where y is the ratio of the energy of the added 
white noise to the energy of the wavelet. 

One can avoid changing R(u)/fTU) at those frequencies where the 
signal-to-nolse ratio Is acceptable by introducing a minimum wavelet 
amplitude or waterlevel (Helmberger and Wiggins, 1971; Clayton and 
Wiggins, 1976), which is expressed as the fraction of the maximum wavelet 
amplitude, |CT(u)| Max ** 
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minimum wavelet amplitude a k]M(oi) ||«| ax (19A) 

where k is called the waterlevel parameter. Expression (19A) 
is used in the denominator of (14) at those frequencies where the 
wavelet amplitudes are low: 

EU) = Wti*) 5 . (19B) 

{max{|W( w )|,k|W( w )| Max }r 

The noise also can be supressed by determining the frequency 
band, (u min » in which the signal-to-noise ratio is acceptable 

and by introducing constant amplitudes outside this band: 



R(u))/W(u) , 

if 

U) . < U) < 

min ~ — 

"Max (20A) 

(a(«) 3 

R Hrin> /ff <W* 

if 

U) < U) • 

min 

(20B) 


R( “Max )/i ' ( “M a x ) ' 

if 

U) < Oi 

Max 

(20C) 


The amplitudes of the calculated input can also be zeroed outside this 
band (Plesinger et al . , 1970): 




R(u))/W(u)) , 


if "min ^ "Max 


if “ < "min or "Max 


(21A) 

(218) 


0 , 


< U) 
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Although there Is not much practical difference between (20) 
and (21) In case of large signal -to-nolse ratio, (20) Is closer to 
reality, since the Input seldom Is a band-limited signal. 

All of these techniques for quality control Involve the 
estimation of the signal -to-nolse ratio, that Is p(w) or Its Inverse 
Fourier transform. This estimation leads to the percentage of white 
noise to be added to the wavelet, the value of y In the y-Fourler 
filter, the waterlevel, or the reliable frequency band, and It requires 
the estimation of n(t) and m(t) or N(u>) and M(o»). The reliable frequency 
band can also be defined by trial and error; the spectral amplitudes 
of the calculated Input would Increase at those frequencies where the 
| l/ff(ui) | becomes small in relation to the noise (Pleslnger et al_., 1970). 

2.3. Instrumental deconvolution of the peaked-mode long-period 
lunar seismograms. 

As shown In the previous section, the noise in the calculated 
ground motion should be suppressed. This noise depends on the noise 
level of the recorded seismograms to be corrected for instrumental 
response and may be estimated from the difference between the observed 
and theoretical calibration pulses. In this section, I describe how the 
calibration pulses were used to estimate the reliable frequency band 
for instrumental deconvolution of the peaked-mode long-period lunar 
sei smograms . 

The simplest and most rapid way to estimate the signal -to-nolse 
ratio of the calibration pulse is by inspection of its phase spectrum 
(see Figure 12 for components 15X and 15Z). Oscillations of the phase 
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spectrum of the observed calibration pulses Increase with decreasing 
slgnal-to-nolse ratio. The phase spectra of the peaked-mode calibra- 
tion pulses are oscillatory at frequencies higher than about 1.3 Hz 
(1.6Hz for component 1 5Z ) (see Figure 12). Near this frequency the 
amplitude of the noise is about that of the signal. The lower limit 
for the reliable frequency band was estimated assuming the additive 
noise to be white, i.e., the signal-to-noise ratios are about the 
same at those frequencies where the signal amplitudes are the same. 
Thus, the lower limit is the frequency where the amplitude response 
of the seismograph is about the same as that at the upper limit and is 
about 0.28 Hz for most components and about 0.17 Hz for component 1 52. 

The estimate for the reliable frequency band for instrumental 
deconvolution was refined by calculating the residuals between the 
theoretical and observed calibration pulses (Kulhanek and Lima, 1970, 
used an analytically calculated amplitude spectrum to estimate digitiz- 
ing errors). The lower curves of Figure 13 show the absolute value of 
the differences between the amplitude spectra of the observed and 
theoretical calibration pulses in the peaked mode. The Fourier 
transforms of the observed and the calculated calibration pulses C"(u) 
and C%), respectively, are given from equations (1), (2), (3), (9), 
and (13): 


C ( co) = c( jcu) ^W(u) + N c (u) 

(22A) 

^(u) = c(jw) ^"(u) 

(22B) 
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where N r U) Is the additive noise In the observed calibration pulse. 

The difference between the Fourier transforms of the calculated and 
observed calibration pulses, aC’U), Is given by equations (22A), (22B), 
and (15): 


* ^(w) - IT(u) » c(j«r 3 M(u) - N (o») (23A) 

V 

And the absolute value of the difference Is 

|«CU)| <LCuf 3 |M(u>)| + |N C (»)| (23B) 

Since M(o>) and N(u>) were assumed to be random, and <^C(o>) was found to be 
random. 


CO) 3 |M(o>)| « |N c (o))| (23C) 

With (23C) , (23B) becomes 

|6C(«)| * |N C («)| (23D) 

Figure 13 also gives the absolute values of the relative errors 
defined as the ratio of the observed-minus -calculated values to the 
theoretical amplitude spectra of the calibration pulse. Table 8 lists 
the upper limit for the frequency band in which the relative errors are 
less than 0.4, 1.0, and 2.5. These relative errors were used to estimate 
the relative difference between the observed and calculated ground 
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motions. If the primary source of noise is the digitizing process and 
the transducer-amplifier system (see Section 1.4.), then the noise is 
signal -independent and 

N(u>) * N (at) . (24A) 

(17D) can be simplified assuming |M(u)|aO and using (24A) and (230): 

| p ( tt )| . Mali (24B) 

|R(u) j 

This equation can be freed from R(u) if g(t) is assumed to be a 
stationary random time sequence. This is ah acceptable assumption at 
those frequencies where the noise amplitudes are comparable to the 
signal amplitudes. This assumption means that 

I s e * T[C(t) [dt ( 25A) 

|R(w)| /jr(t)Jdt 

(25A) gives the simplest error estimation using (24B); 

|p(u) I » iMmlLc (25B) 

|Cl«)l 

The least-squares technique of fitting the calibration pulse 
with the theoretically calculated pulse also yielded relative noise 
parameters, such as the maximum and the average of the residuals in the 
frequency band of 0 to divided by the amplitude spectrum of the 
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calibration pulse at the peak frequency. The schematic flow chart for 
calculating these values Is shown In Figure 17. These relative noise 
parameters can be used as waterlevels If they were corrected for the 
amplitudes and the noise of the records which may be different from those 
of the calibration pulse (compare with equation (25B) ) . Since most of 
the long-period lunar seismograms are weaker and noisier than the 
calibration pulse for component 15Z and are more energetic than the 
calibration pulses for the other components, two waterlevel parameters 
were chosen: 0.03 for seismograms with amplitudes larger than 30 to 50 
digital units and not too noisy and 0.1 for the smaller signals. The 
reliability of these waterlevels was tested by deconvolving some- of the 
records with different waterlevel parameters before the routine process- 
ing. These trial deconvolutions showed that the above choices were 
satisfactory except for the very noisy records. 

The noises which are easily spotted and classified by inspection 
are the following: 

1. Most of the "glitches" are those Impulsive- type noises 
In the peaked mode which look like distorted calibration pulses and 
can be described rather accurately by the following empirical formula 

CjC(t) + c 2 ^I * c//{/c ] w(t)dt + c 2 w(t)}(dt) 2 (26) 

dt 

where c^ and c 2 are constants, C(t) is the theoretical or observed 
calibration pulse, C Is given by equation (9), and w(t) is the impulse- 
displacement response of the seismograph. The times of occurence of 
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these glitches correlate strongly with temperature changes near the 
station. A glitch was easily removable by the model described by 
equation (26) If the power spectrum of the seismic signal was stronger 
at frequencies above the peak frequency of the seismometer than at the 
peak frequency. However, If there were many overlapping glitches and 
the seismic signal was strong near the peak frequency of the seismograph 
response, removing the glitches by fitting the seismograms with equation 
(26) distorted the signals considerably. This distortion was caused by 
not having been able to resolve the exact beginnings of the individual 
glitches. 

Other techniques, such as the prediction error operation (Wiggins 
and Miller, 1972) and adaptive filters (Sims and D'Mello, 1978), were 
not tried In this study. They also would distort the signal, because 
the distortion increases rapidly with the unpredictability of the 
glitches in comparison to the unpredictability of the signal (Wiggins 
and Miller, 1972). 

2. Bit-holds and bit-jumps of the analog-to-digital converter 
were observed only at station 14. These types of digitizing errors are 
Illustrated in Figure 22. It is clear from Figure 22 that the bit-holds 
cannot be removed from the seismograms and that the bit-jumps can be 
removed without any information loss. The correction for bit-jumps was 
done by checking six-minute sections of the seismograms for missing 
levels of digitization. 

3. Most of the spikes were also removed from the seismograms 
by checking the changes in signal amplitudes in. windows of few seconds. 
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The Instrumental deconvolution was carried out In the frequency 
domain as the schematic flow chart of the computer program shows In 
Figure 19. After a time section of the seismogram was corrected for 
bit-jumps and for spikes, the d.c. level was removed, and cosine-bell 
tapers were applied at both ends of the time section. The Fourier 
transform of this tapered signal was divided by the peaked-mode complex 
transfer function of the seismograph according to (19). The transfer 
functions were adjusted with the waterlevels by a separate computer 
program whose schematic flow chart Is shown In Figure 18. The 
Instrumental deconvolution was completed by taking the Inverse Fourier 
transform. The amplitude spectra of the deconvolved records reveal 
that other important filtering effects characteristic of each recording 
site remained in the lunar seismograms. These filters cause most of 
the amplitude spectra to deviate systematically from a rather uniform 
spectrum and prevent an easier comparison of the different seismograms. 
The removal of these effects is discussed in the following section. 
Naturally the approach cannot be deterministic, as correcting the 
seismograms for the instrumental response, and should rely on statistical 
means. 


2.4. Design of the whitening filters to remove the filtering effects 
of the near-surface zone. 

In addition to the effects of the peaked-mode Instrumental 
response, the following major processes also Influence the observed 
amplitude spectra of the long-period lunar seismograms. 


1. Source effects 



2. Propagation effects, e.g., frequency-dependent attentuatlor. 

3. Noises, e.g. spikes and glitches 

4. Near-surface effects, e.g., large velocity gradients, 
lateral Inhomogeneities 

5. Coupling between the seismometer and the ground. 

Although all of these effects are Important In shaping the lunar 
seismograms, the filtering effects of the near-surface zone and the poor 
coupling are considered to be primarily responsible for most of the 
differences observed among the signals recorded by the various components 
(lamnleln et aL , 1974; Nakamura et aK , 1975; Mark and Sutton, 1975). 
Apart from scaling the signal amplitudes of all three long-period 
components of the ground motion to the same level (Goins, 1978) the 
seismograms were not corrected In the past for these effects. 

Since no deterministic way is available to reduce these filtering 
effects one must resort to statistical methods. Ideally, the statistical 
approach would mean choosing a set of events for each similar group of 
sources such that the hypocenters are evenly distributed In space and the 
noise level Is rather low. It would also mean combining the amplitude 
spectra properly within each group; that Is, the statistical combination 
of the indlvludal amplitude spectra would result In cancellation of all 
Influences except those attributed to the near surface and the source 
effects . 

Initially, I selected 65 of the largest events from the complete 
seismic data set. These 65 events Include 40 Impacts (8 artificial and 
32 natural), 10 shallow moonquakes, and 15 deep moonquakes. At least 
one time section of 2048 samples (309 seconds) was chosen from each 
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record starting just before the beginning of the seismic signal. For 
signals from some distant events, another section of the same length, 
following the first by 3 or 4 minutes, was also selected. These 65 
events and the beginnings of the approximately five-minute sections are 
listed In Table 9. The first letter of the code In this table refers 
to the type of source: A to deep moonquakes, H to shallow mocnquakes 
(also called high frequency teleselsmlc events), C to Impacts. The two 
digit numbers following the first letter are used to denote source loca- 
tions for deep moonquakes and are arbitrary sequence numbers for shallow 
moonquakes and Impacts; A46 was reclassified as A01 by Nakamura (1978). 
The artificial Impact records are distinguished from the natural Impact 
records by using a letter and a number following the first letter. The 
letter denotes the source (L for Lunar Module and S for Saturn IV 
booster), the number Is the second digit of the Apollo mission number 
(12, 14, etc.). 'WL' Is the waterlevel parameter In percent used In 
the peaked-mode Instrumental deconvolution process (see Section 2.3.). 
'Com' gives information about the various components for each station, 
such that Integers 1, 2, or 3 imply that the X, Y, or Z component, 
respectively, did not operate properly for that event and thus no signal 
can be seen; the letters X, Y, or Z are used to denote those time 
sections whose amplitude spectra peaked In the frequency band of 
approximately 0.3-2 Hz. Since this selection eliminated those time 
sections which possess Inadequate signal -to-noise ratios, sections 
marked with X, Y, and Z were used for designing the whitening filters. 

Table 10 lists the number of time sections used In this study 
by components and by event types. Those listed In Table 10 as 'Good' 
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met the above slgnal-to-nolse ratio criterion; those In column 'Noisy' 

Mere not used to design the whitening filters. Table 10 also lists the 
number of time section? which were recorded In the flat mode. 

To obtain the relative magnifications for the various components 
at each station , the same events should be chosen for each component. 

Since the small number of time sections meeting the above slgnal-to-nolse 
ratio criterion does not permit this, It was necessary to compromise by 
selecting the same events for the horizontal components only, except for 
the few events listed In Table 9. The present selection of events per- 
mits rotation of the horizontal components Into radial and transverse 
directions and particle-motion studies. 

It was also likely that unbiased statistical analysis cannot be 
performed with the small number of events within each of these event 
types. It was reasonable, nevertheless, to expect that the seismic 
records can be corrected for the average near-surface effects at each 
site If the other effects, such as the spectral density functions of the 
Individual sources and the propagation effects along the various raypaths, 
produced smaller alterations In tie spectral amplitudes than the near- 
surface effects. To test this supposition, amplitude spectra of time 
sections from the three major source types were combined at each site as 
weighted arithmetical averages. 


T s (nc) 


’sH 

Z c,X„.,(na) 


n sh i:r™ 


(27A) 
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7 (no) ■ J- r c.Y >(no) 
s n sH 1*1 1 S1 


N 


sH 


(27B) 


7 s (no) 


N sZ 

u 2 CiZ i(flo) 
N sZ 1-1 1 81 


(27C) 


where n ■ 0,1 ,2... . ,1023; superscript s denotes the source types, 

X s1* Y s1, and Z $ ^ are the amplitude spectra of the Individual time 
sections, recorded In the peaked mode and corrected for the Instrumental 
response, X x , 7 , and I s are the average of the and N x2 number of 

Individual amplitude spectra, o Is the sampling Interva' In the angular 
frequency domain (o - 2ir* 0.003235 sec”^), and are individual weight 
factors . 

The flow chart of the computer program for calculating 
spectral amplitude averages according to equations (27) Is shown In 
Figure 20. The averaging can be done after the Individual amplitude 
spectra have been smoothed and, if desired, weighted. The smoothing 
will not be considered at this point since the associative law holds for 
summing and smoothing with moving averages In rectangular windows. When 
no weighting ic applied to the Individual spectral amplitudes (c^l), 
the seismic records with Targer amplitudes have greater weight In the 
average amplitude spectrum than the weaker records. When the energies 
or the amplitudes of the sections are scaled to the same level before 
simulation, the Individual amplitude spectra are treated with approxi- 
mately the same weight, except that a noisier record contribute less in 
the frequency band of Interest than the records with good signal -to-nolse 
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ratio If the noise is stronger outside this band. Weighting does not 
permit derivation of relative magnification ratios for the two hori- 
zontal components at a particular site, but allows an empirical 
reliability test for the consistency of the shapes of the determined 
average amplitude spectra. 

In addition to this test for the consistency of the average 
amplitude spectra, subsets were selected within the most numerous 
groups. The average amplitude spectra of these different groups are 
similar to each other for- a given component except In certain frequency 
bands in which noise, discussed in the following section, is obviously 
present, and except for impact signals which are most affected by distance 
dependent attentuation. Based upon the similarity among the average 
amplitude spectra of the different groups and upon the similarity between 
the amplitude spectra of two consecutive five-minute sections, I 
conclude that the five minute sections starting at about the beginning 
of the seismic signals adequately defines the statistical properties of 
the amplitude spectra. The average amplitude spectra obtained by 
simple addition of the amplitude spectra for the different components 
and source types are shown in Figures 23 through 28. 

Although there are some differences among the average amplitude 
spectra of the different source types, I stacked these spectra for every 
component because the number of sections is small for some of these 
groups and I wanted to derive only one whitening filter for each compon- 
ent. These average amplitude spectra are shown In Figures 29 through 32 
after they were smoothed with moving averages In a rectangular window of 
17 samples (0.055 Hz). Figures 29 through 32 also show the average 
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horizontal amplitude spectra at each station defined as 

HOR(na) *\^T(no) 2 ♦ 7(no)\ n»0,l ,2,. ...1023 (28) 

where J and T are the average amplitude spectra for components X and Y, 
respectively. 

These average amplitude spectra can be considered as the 
amplitude spectra of the near-surface effects for an average source 
and attentuatlon. The whitening filters can be obtained from the 
smoothed average amplitude spectra assuming the filters to be either 
zero phase or minimum phase. In the latter case, the filter is causal 
and the phase spectrum can be calculated with the Hilbert transform of 
the logarithm of the amplitude spectrum. The minimum-phase assumption 
contradicts the fact that the signal was dominated by scattered surface 
waves and the zero-phase assumption violates causality. However, I 
adopted the assumption of zero phase in calculating the whitening 
filters because the true phase spectra of the direct shear and secondary 
arrivals were not expected to be determined from the amplitude spectra 
of the long time-sections (five-minute sections were required to 
describe the average spectral properties of the scattered surface waves). 
To avoid the undesirable sidelobes of a rectangular window, the average 
amplitude spectra were smoothed with moving averages over a 0.109 Hz 
window which was cosine tapered at both ends around the central fre- 
quency of the window. Thus, the amplitude spectra of the whitening 
filters differ slightly from those of Figures 29 through 32. 

To reduce the numerical errors, the zero phase whitening filters 
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were multiplied by the complex transfer function of the seismograph and 
equation (20) was used to stabilize the Inverse filtering process In 
the peaked mode. (Figure 21 shows the flow chart of the computer 
program). The reliable frequency band was chosen to be (0.2 Hz, 1.5 
Hz). The upper limit of this band could have been slightly higher 
for a few of the largest signals. Seismic signals recorded In the flat 
mode were also corrected with the whitening filters. This Inverse 
filtering process did not require suppression of the noise because the 
whitaning filters contained some noise thought to be white. The 
whitening filters in the peaked mode could have been designed such 
they also Included the instrumental effects. Although this approach 
would have eliminated the necessity of determining the transfer function 
of the seismographs, it would not have yielded information of the 
seismographs, of the near-surface structure, and of the noise, and 
would not have permitted correction of the flat-mode seismograms for the 
near-surface effects. 

Although the zero-phase nature of the whitening filters caused 
the seismic arrivals to be observed slightly earlier, the Inverse 
filters diminished the sharp peaks In the amplitude spectra considerably. 
Tests conducted In the time domain also showed that the inverse filters 
made the comparison among the different components possible. This is 
demonstrated In Figure 33 with the deconvolved Y components of ground mo- 
tion at stations 14 and 16. The knowledge of the average amplitude spectra 
also permits derivation of the average relative amplification between the 
horizontal components, which is discussed in the following section. 
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2.5. Some properties of the amplitude spectra of the long-period 
lunar seismograms. 

The average amplitude spectra derived In the preceding sections 
permitted empirical determination of the relative amplitude ratios between 
the horizontal components of the long-period lunar seismographs. 

The ratios of the average horizontal amplitude spectra (X/Y) 
of all sections with good slgnal-to-nolse ratio are shown in Figure 34 
after they were smoothed with moving averages in a rectangular window 
of 17 samples (0.055 Hz). The relative average amplifications calcu- 
lated for two frequency bands are listed In Table 11. Although the 
X/Y ratios at some frequencies differ from the average by as much as 
30% and might indicate the presence of lateral inhomogeneities, the 
differences among the amplifications at the different stations are 
primarily attributed to the coupling between the seismometer and ground. 
The cable connecting the sensor unit to the central electronics unit 
is thought to be the most important factor in the coupling (Lammlein 
et al_. , 1974). This cable oomes out at an angle of about 30° to the 
positive Y direction and possibly acts like a fourth leg. Sharp 
increases in the amplitude spectra at the peak frequency of the ampli- 
tude responses of the horizontal components at station 16 are also 
interpreted to be a result of this mechanical effect. 

At certain frequencies most of the time sections show unusually 
high amplitudes which are associated with noise. The large amplitudes 
at less than 0.2 Hz correspond to the frequency content of the 
glitches as mentioned in Section 2.3. They can be buried in the signal 
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and unnoticed in the seismograms recorded in the peaked mode. These 
glitches are abrupt ground motions caused by thermal expansion of the 
temperature insulating shroud and of the soil surrounding the 
seismometer. 

Other types of instrumental noises, such as bit-jump and bit-hold 
errors of the analog-to-digital converter and spikes, were discussed in 
Section 2.3. Some of the long-period signals show a spectral peak 
at the frequency of 1.6563 Hz or half of the Nyquist frequency, which 
corresponds to one data frame (see Section 1.2. for details about 
instrumentation) and is probably caused by varying reference voltage 
in the analog-to-digital converter. Other noises of possibly instru- 
mental origin are present in the seismic records from stations 12 and 14 
at periods of 3-4 seconds. Station 12 seismograms also show a spec- 
tral peak at periods of 8-10 seconds. Both of these noises are 
attributed to the on/off switching of the temperature control unit. 

The mechanism of this interaction is unknown. 

2.6. Summary of Part 2. 

This part of the dissertation has described the methods used for 
correcting the long-period lunar seismograms for instrumental response. 
To control the noise, water! evel parameters of 0.10 for most seismo- 
grams and of 0.03 for seismograms with good signal-to-noise ratio were 
determined from the analysis of the signal-to-noise ratio of the 
calibration pulses. The instrument-corrected seismic records were used 
to determine the whitening filters for each component in order to 
reduce the characteristic near-surface effects and to make uniform 
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filtering of the different seismic records possible. These filters, 
which can be considered to be the amplitude spectra for an average 
source and propagation character, work satisfactorily. Although the 
zero-phase nature of the whitening filters causes the seismic arrivals 
to be observed slightly earlier, the sharp spectral peaks of the long- 
period lunar seismograms are diminished considerably. The rather white 
nature of the amplitude spectra permits comparison of the seismograms 
from different stations as described In Part 4. The average spectral 
amplitudes were also used to derive relative amplitude ratios between 
the horizontal components. , 



PART 3. NEAR-SURFACE SHEAR WAVE VELOCITIES AS INFERRED FROM THE 
SPECTRAL AMPLITUDE RATIOS OF THE LONG-PERIOD LUNAR SEISMOGRAMS. 

3.1. Introduction. 

In the preceding part of this dissertation, the amplitude 
spectra of the most energetic long-period seismograms were studied In 
309-second sections. Although the amplitude spectra depend on the 
source depth and distance, they are similar In two consecutive sections 
of the signals. The average amplitude spectra calculated for the three 
major types of sources at a given site also show a high degree of 
consistency. These observations suggest that the average amplitude 
spectra may be interpreted in terms of near-surface effects. This 
interpretation is important because other information on the physical 
properities of the near-surface zone is very limited. 

Similar congressional velocities of 92 to 114 m/sec were 
obtained for the topmost layer, the regolith, by active seismic experi- 
ments at the Apollo 14, 16, and 17 sites (Cooper et ajL , 1974) and by 
analyzing the signals from the Lunar Module lift-off at the Apollo 11, 
12, 14, and 15 sites (Nakamura et al_. , 1975). Cooper et al_. (1974) 
characterized the regolith as a rather homogeneous layer of differing 
thicknesses at the different sites (8.5, 12.2, and 4.0 m at stations 
14, 16, and 17 respectively) and determined the compressional velocities 
to be about 300 m/sec for the underlying material at all three sites. 

The results of the seismic profiling experiment at the Apollo 17 site 
gave velocities of 495, 960, and 4700 m/sec at depths of 32, 390, and 
1385 m (Cooper et al_. , 1974). The last two velocities are less certain 
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because the few arrivals measured at distances greater than l.i km 
had to be corrected for topography. Although the velocities were given 
as stepwise Increases with depth, the Investigators recognized the 
possibility of a continuous Increase of the velocities with depth. In 
fact, the travel times and amplitudes at the Apollo 14 and 16 sites 
were re- Interpreted In terms of a power-law velocity variation with 
depth, v(z) » 110z^® m/sec, where z Is the depth In meters down to 
10 meters (Gangl, 1978). 

These low velocities and their rapid Increase with depth were 
shown by laboratory measurements to be reasonable for porous, 
brecclated material with no volatile content. Measurements by Stesky 
and Renton (1977) on lunar soil samples under low loads yielded 
velocities appropriate for the uppermost 100 m and suggest that the 
compresslonal velocities and the Poisson's ratio did not depend 
strongly on the density of the soil samples. The Poisson's ratio is 
about 0.46 at near-zero pressure and Is rather constant at 0.40-0.43 
In the pressure range of 0.2 to 5 bars. The compresslonal velocity 
increases from 250 m/sec at the near-zero pressure to 750 m/sec at 
5 bars. Laboratory measurements on Apollo 17 rock powders by Talwani 
et al_. (1973 and 1974) yielded compresslonal velocity gradients between 
0.5 and 0.1 (km/sec)/km at confining pressures of 0 to 0.5 kbars, and 
10 to 20 times less at higher pressures of up to 2.5 kbars. The 
velocities and their gradients are higher for the "light mantle 
material", which Is a product of a variety of breccias, than those for 
the "dark mantle material", which has a composition similar to that of 
mare basalt samples (Talwani et al., 1974). 
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Soil densities were measured in the laboratory (Carrier et al_. , 
1973; Stesky and Renton, 1977) and by soil mechanics experiments in situ 
(e.g., Mitchell et ajL, 1972; Houston et al,. , 1974). The results, which 
were summarized by Helken (1975), show widely varying densities of 
less than 1 to more than 2 g/cm . Sutton and Duennebler (1970) obtained 
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densities of 1.1 to 1.5 g/cm at the landing sites of Surveyor space- 
craft by modeling the influence of the soil on spacecraft vibrations. 

The ratios of the horizontal to vertical long-period spectral 
amplitudes were studied by Nakamura etal,. (1975) and by Mark and Sutton 
(1975). The former group obtained the thickness and the shear velocity 
of the regolith from the maximum of the spectral amplitude ratios, 
assuming that it was produced by quarter-wavelength interference of 
shear waves trapped in the regolith. An average shear velocity of 
37 m/sec was found in a layer of 4.4, 3.7, and 4.4 meters thick at the 
Apollo 11, 12, and 15 landing sites, respectively. Mark and Sutton 
(1975) fitted the spectral ratios of a few impact records with spectral 
ratios calculated theoretically for the fundamental -mode Rayleigh waves. 
They found that the shear velocities- increase from about 35 m/sec at 
the surface to about 400 m/sec at a depth of 120 m under stations 12 and 
15. It was also shown that the velocity stratification obtained at the 
Apollo 14 site by Cooper et al_. (1974) results in a theoretical infinite 
value near the peak of the observed spectral ratio curve. 

In this part of the dissertation I also used the spectral ratio 
technique, modeling the observed amplitude ratios with ellipticities 
calculated for the fundamental -mode Rayleigh waves, to obtain the shear 
velocity distributions at the Apollo 12, 14, 15, and 16 landing sites. 
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First, the spectral amplitude ratios for many of the largest events are 
shown to be similar at a given site. Then, the average spectral 
amplitude ratio curves are fitted with theoretical curves. A discussion 
of the velocity distributions thus obtained concludes this part of the 
dissertation. 

3.2. The observed spectral amplitude ratios. 

The spectral amplitude ratios, IT, were calculated from the 
average spectral amplitudes (see equations (27) and (28)) according to 
the following relation 


1T( ro ) . HOR(no) . Vjjdl t T(na]j 
T(na) Z(na) 


n=0,l,2 1023 (29) 


where Y, 7, and 7 are the average amplitude spectra calculated from the 
spectral amplitudes of all time sections having adequate slgnal-to-nolse 
ratios. These spectral -ratios are shown as the top curve of Figures 35 
through 38 after they were smoothed with moving averages in a rectangu- 
lar window of 17 frequency samples. If the attenuation and source 
effects were the same for the vertical and horizontal components and 
the same events were used, then these average spectral ratios would be 

free of such effects. Since the events used In calculating the average 
horizontal and vertical amplitude spectra (see Tables 9 and 10) are 
different, the spectral ratios calculated by (29) may Include effects 
other than those of near-surface origin and may be distorted at least by 
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a scale factor. To test the Importance of these effects on the spectral 
ratios and to obtain the correct values for these ratios, the spectral 
ratios were calculated according to the following equation for those 
time sections which met the slgnal-to-nolse ratio criterion described 
In Section 2.4. 

Vw ^) 2 + Y >°) 2 

R .(no)-— - , n-0, 1,2,. ...1023 (30) 

Zj(no) 

where R^ are the Individual spectral ratios, XJ, Yj, and Z} are the 
spectral amplitudes In the horizontal and vertical directions after 
they were smoothed with Hammlng-Tukey smoothing coefficients (BSth, 

1974, p. 179). To obtain a better estimate for the spectral ratios, 
these individual spectral ratios were averaged for each source type: 


R s (no) 


E R^ (na) , 
N s 1-1 1 


n-0,1 ,2, ...,1023 


(31) 


where subscript s stands for the source type and N s is the number of 
sections. The flow chart for calculating these average spectral ratios 
Is given in Figure 20, and the average spectral ratios are shown as the 
lower three curves of Figures 35 through 38. 

As seen from these figures, the average spectral ratio curves 
calculated by (29) and (31) are very similar at a given station, except 
those at station 15 and that for deep moonquakes at station 12. At 
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station 15, the dominant peak In the average spectral amplitude ratio 
curve Is apparently shifted toward lower frequencies with Increasing 
source depth, while at station 12, there Is no dominant peak In the 
spectral ratio curve for deep moonquakes. The moonquake signals are 
generally sraller than those of Impacts, and stations 12 and 15 have 
smaller horizontal amplitudes than stations 14 and 16. Since the 
spectral amplitude ratios are dominated by the noise at those frequen- 
cies where the signal amplitudes and/or the Instrumental sensitivities 
are low In comparison with the noise level, the apparent source depend- 
ency of the spectral ratios may be attributed to the use of less 
energetic signals recorded by the very sharply peaked seismometers at 
the least sensitive stations. Indeed, the spectral ratio curves of 
the largest moonquakes are much more like those of Impacts than those 
of weak moonquakes . 

Thus, the spectral ratios at a given station are generally 
independent of the source depth and distance. Furthermore, the spectral 
ratios of the two different sections of the signals from the distant 
events are quite similar. The hypothesis that spectral ratios are 
determined primarily by the physical parameters of the near-surface 
zone fits these facts very well. Scattering and ground coupling may 
also influence the spectra ratio curve. 

The problem of comparing the measured spectral ratios for the 
Isotropically arriving scattered surface waves with theoretical ratios 
for plane Rayleigh waves was solved by Introducing a multiplying factor 
of 2»^/tt»0.9 to equation (30) (Mark and Sutton, 1975). However, I 
omitted this factor because the relative amplifications among the 
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components are not precisely known (see Part 1 of this dissertation) 
and Imperfect coupling between the seismometer and the ground is si so 
suspected of causing certain Irregularities, such as the sharp Increase 
of the station 16 spectral amplitudes at 0.45 Hz (see Section 2.5.). 
Since the signal amplitudes of the Y component of ground motion are 
about 1.7 times larger than those of the X component for station 16, 

(see Section 2.5.), and since the sharp peak of the spectral ratio 
curve at about 0.45 Hz corresponds to the peak magnification of the Y 
component operated in the peaked mode, the observed spectral ratios at 
around this frequency were excluded from the frequency band in which 
the observations should be fitted with theoretical curves. The lower 
limit for this frequency band Is the frequency at which the long-period 
noise is sufficiently small and is about 0.3 Hz for all stations. The 
higher limit is the frequency where the horizontal -to-vertical spectral 
ratio becomes that of the noise, i.e., 2.3, 1.6, 2.1, and 1.7 Hz for 
stations 12, 14, 15, and 16, respectively. 

3.3. The theoretical spectral amplitude ratios calculated for the 
fundamental -mode Rayleigh waves. 

Assuming the vertical velocity variations to be much larger than 
the lateral ones, the near-surface structure was modeled as a series of 
homogeneous, isotropic, horizontal layers. The validity of this model 
can be partially tested by calculating the ratio of the average X and Y 
amplitude spectra for all sections having good signal-to-noise ratio. 

The X/Y spectral amplitude ratios were discussed in Section 2.5. and are 
shewn in Figure 34. Some of these ratios differ from the average 
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(listed In Table 11 for two frequency bands) by as much as 30%, and 
follow a pattern similar to the horizontal -to- vertical spectral 
amplitude ratios near the dominant peak of the horlzontal-to-vertlcal 
spectral ratios. These variations may Indicate a tilt In the particle 
motion ellipse, which may be produced by anelastlc effects of the 
medium (Boore and Toksoz, 1969), by the presence of anisotropy, by' 
lateral Inhomogenltles, or by ground coupling effects (e.g., a rocking 
motion of the Instrument). Variations of the X/Y spectral ratios are 
not considered to be dominant effects. 

Assuming a horizontally layered medium, the ratio of the 
horizontal and. vertical axes of the particle motion ellipse, also 
called elllptlclty of the Rayleigh type surface waves, was calculated 
theoretically with the Haskell matrix technique (Haskell, 1953). Since 
the distribution of energy among the different Rayleigh modes Is 
unknown, I also assumed that Rayleigh waves of the fundamental mode 
dominates the ground motion In the frequency band of Interest. Disre- 
garding the possible contribution of body waves, Love waves and higher- 
mode Rayleigh waves appears to be justified quantitatively by the good 
agreement between the results presented here and those from refraction 
experiments, as described In the following section. Some additional 
qualitative arguments can also be given. The long time-windows used and 
the similarity between the spectral ratios In two consecutive five-min- 
ute windows support the assumption that the contribution of the body 
waves is negligible. If the energy of Love waves varied with frequency 
In the same way as the energy of the horizontal ground motion of Rayleigh 
waves, then the observed spectral ratios would differ from those 
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calculated for Rayleigh waves by a scale factor only. The low veloci- 
ties and large velocity gradients near the surface, as well as 
heterogeneities concentrated near the surface, support the assumption 
that higher-mode Rayleigh waves may be neglected because surface sources 
produce fundamental -mode surface waves primarily and because the 
conversion from fundamental -mode surface waves to higher-mode waves is 
small in a scattering medium (Mai in, 1978). To emphasize the contri- 
bution of the fundamental -mode Rayleigh waves and to use the most 
energetic events, only signals from meteroid impacts were used. 

The densities and the velocities should satisfy certain con- 
straints established by other investigations, which were reviewed in 
Section 3.1. In order to simplify the comparison among the different 
sites and to handle as few variables as possible, I varied the thickness 
and the sheai velocity of the layers and kept the compressional-to-shear 
velocity ratios and densities identical in the corresponding layers at 
all sites. These restrictions are also justified by the fact that the 
theoretical ell ipti cities of the fundamental -mode Rayleigh waves in the 
0.3-2. 3 Hz band are most sensitive to shear velocities and least to 
compressicnal velocities and densities. Since the theoretically calcu- 
lated el 1 i pti ci ties are invariant to a corrcno" scale factor simultaneously 
applied to the thicknesses and velocities (Haskell, 1953), the shear or 
compressionai velocity must be fixed at a given depth in order to obtain 
unique solutions. A compressiona. velocity of 104 m/sec in the 
regolith, v.f.ich is the a''t.age of the observations at the four sites, 
was used to fix this scale factor. Moreover, the density of the super- 
ficial layer also was fixed because the theoretically calculated 
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elllptlclties are affected by the density ratios only. Finally, It 
was also my objective that the theoretical models consisted of as few 
layers as were necessary to model the spectral amplitude ratios. 

Two different compresslonal-to-shear velocity ratio and density 
profiles were assumed for calculating the theoretical elllpticltles. 

In the first model, model 'A', which Is similar to the model of Mark 
and Sutton (1975), the compresslonal-to-shear velocity ratios are 
higher and the densities are lower in the four uppermost layers than In 
model * B ' . Average spectral ratio curves for Impacts were fitted with 
theoretical curves by adjusting the thicknesses and the velocities of 
the layers of the model. These adjustments were made by examining the 
partial derivatives of the theoretical el 1 iptlci ties with respect to 
these variables (a least-squares iterative scheme was also tried, but 
the Infinities in the calculated ellipticitles caused the Iteration to 
diverge). The fit in the frequency band of Interest was judged by 
eye on plots having linear frequency and amplitude scales, such as 
those of Figures 35 through 38. The spectal ratios calculated for the 
final models are shown In Figures 35 through 38 as a dot-dash line for 
model 'A' and a dashed line for model ' B ' . 

The structural parameters of these models are given in Tables 12 
and 13. The agreement between the theoretical and observed ratio curves 
is reasonably good except at high frequencies, where noise dominates the 
spectra, and at the low frequency end of the spectral ratio curves for 
stations 14 and 16. The sharp increase in the station 16 spectral 
ratios at 0.45 Hz is attributed to a rocking motion of t» seismometer, 
as was discussed earlier. The decrease below 0.8 Hz at station 14 
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could not be modeled without the introduction of Infinities in the 
theoretical curves. 

3.4 The near-surface velocity distributions. 

The discontinuities introduced to facilitate the theoretical 
calculations do not necessarily represent actual interfaces. Although 
I cannot rule out the presence of discrete layers, such as the regolith, 
the velocities of models 'A' and ' B * , which are remarkably similar, may 
be interpreted as varying conti nous ly with depth. Figure 39 shows the 
shear velocity distributions of model * B ' . This model, which possesses 
smaller compressional-to-shear velocity ratios and higher densities, is 
preferred to model 'A' because its theoretical spectral ratio curves fit 
the observed curves more closely, the velocities at the various sites 
differ less from one another, and the theoretical arrival times of the 
compressional waves (Figures 40 through 44) are closer to the first 
arrivals from the Lunar Module lift-offs (Nakamura et al_. , 1975) at 
stations 12, 14, and 15. 

All model '8' velocity distributions are similar and differ 
from the average for all sites by not more than 20%, except for the 
bottom layer at station 14. The velocities at stations 14 and 15 are 
almost the same down to about 50 meters, below which the velocities are 
much higher at site 14 than at 15. The thicknesses of the layers having 
the same velocities are about 20% larger at station 15 than at 12 in the 
uppermost 41 meters, below which the shear velocities are about the 
same at these sites. The velocities at station 16, except for the 
topmost layer, which is fixed, are lower than those at any ot'.or station. 
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If the calculated velocities are taken as the averages in ar\y 
layer and are assigned to the center of this layer, marked with crosses 
in Figure 38, then a piecewise linear increase of the shear velocities 
with depth in two segments, as shown with dashed line in Figure 38, 
describes the shear velocities of model ' B ' rather well. Each segment 
of these continuous shear velocities, s, is described as 

3 * a + bz (32) 

where z is depth in meters, and a and b are constants, which were 
determined by least squares and are listed in Table 14. Table 14 
also lists the depths where the gradient, b in equation (32) changes, 
i.e., where a second-order discontinuity occurs. Although these depths 
and a in (32) are affected by the same scale factor which was used to 
fix the compressional velocity of the regolith, the value of b in (32) 
is not affected by the choice of the surface velocity. 

The shear velocities obtained by Mark and Sutton (1975) are 
higher in the uppermost 80 to 100 meters than those of models 'A' and 
’ B ' at station 12. The differences are attributed to the way the fit 
was obtained. The details of the spectral ratio curve at short periods, 
which are controlled by the physical parameters of the shallower layers, 
weigh little if fitted on plots of the logarithm of the spectral 
amplitude ratios versus period, as was done by Mark and Sutton (1975). 
They specified the velocities and the density at site 12 down to 413.5 
meters, which was not necessary in this study. Their velocities are 
closer to those of models 'A' and ' B ' at station 15 than at 12. 

Although power-law functions were found to fit the shear 
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velocities of model ' B ' rather poorly, the congressional velocity, a . 
can be described adequately by a power- law function in the 
uppermost 10 to 13 meters: 

a * cz p (33) 

where z is depth in meters, and c and p are constants, which were also 
determined by least squares and are listed in Table 14. The exponent 
in (33) is also invariant to the scale factor applied to fix the 
compressional velocities at the surface. Table 14 also lists the 
number of layers, where this power-law fit is adequate, and the average 
compressional velocities in these layers. 

Compressional velocities of model ' B 1 agree very well with those 
obtained from the seismic signals of the Lunar Module lift-offs 
(Nakamura et al_. , 1975). The first arrivals at stations 14 and 15 are 
within 10 msec of the theoretical arrival times of refracted compress- 
ional waves calculated for model ' 8 ' (Figures 41 and 42) and within 
80 msec at station 12 (Figure 40). 

The average compressional velocity calculated with Gangi's 
(1978) power-law formula is about 134 m/sec in the uppermost 10 m 
which is close to those found in this study at stations 14 and 16 
(see Table 14). From equation (33) the travel times, t, of the com- 
pressional wave arrivals at distance x is 

t * ^(x/Xq) 1 "* 5 (34A) 

where t is the travel time at distance x„ (Gangi, 1978). Obtaining p 
o 0 
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from the best fitting power-law models, p*l/3 (Table 14), and using 
the first arrival from the Lunar Module lift-off at station 14, 
t Q »670 msec and x Q *178 meters (Nakamura et al., 1975), 

t * 21.2’x^ 3 msec (34B) 

where x is in meters. Table 15 compares the travel times calculated 
for various velocity models with those observed from the Apollo 14 and 
16 stacked data. As seen in Table 15, travel times calculated with 
equation (34B) agree best with the observed ones. 

The compressional velocities of models 'A' and ' B ' increase 
more smoothly with depth than those obtained by Cooper et al_. (1975) 
at the Apollo 14, 16, and 17 sites. They obtained a compressional 
velocity of 4.7 km/sec at a depth of 1385 m at site 17. If the 
gradient of the shear velocity for the bottom zone (see Table 14) 
is extended to a depth of a few kilometers and the compressional -to- 
shear velocity ratio is assumed to be the same (2.0), then the extra- 
polated compressional velocities are much higher than theirs between 
depths of 0.1 and 1.0 km and reach 4.7 km/sec at depths of 1250 to 
1500 meters at stations 12, 15, and 16. The differences at depths 
less than 1 km are attributed to the large uncertainties in the 
travel time data of the seismic profiling experiments at distances of 
1.1 to 3 km. The linear extrapolation of the compressional velocity 
of model ' B ' at station 14 yields a value of 4.7 km/sec at a depth of 
790 m. These comparisons suggest the velocity gradients of Table 14 
may be correct in the uppermost one or two kilometers at all sites 
except at 14, where the velocity gradient obtained in this study 
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should decrease with depth. 

The unusually high velocities at depths between 50 and 100 
meters at station 14 also can be supported by the following argument. 

If the spectral peaks observed at 0.85 Hz for the horizontal component 
of the ground motion and at 1.6 to 1.7 Hz for the vertical component 
(see Figure 30) were caused by resonance of the shear and congressional 
waves In a layer, respectively, then the peak frequencies must be 
proportional to the velocities In this layer and, therefore the layer 
must have an average compressional-to-shear velocity ratio of about 2.0. 
Since the compressional-to-shear velocity ratio was found to be much 
higher than 2.0 In the regolith (see Section 3.1), the spectral peaks of 
the ground motion at station 14 must be explained by a deep velocity 
discontinuity. 

3.5 Discussion and conclusion. 

Results from this study indicate that processes acting at the 
surface of the moon have produced a surface zone of remarkably uniform 
physical properties. Although the presence of discrete layering is 
not ruled out by these results, the obtained layered velocity-depth 
functions can be well fitted with two piecewise linear segments 
Continuous velocities, thus derived, are characterized by two zones; 
a thin (10 to 25 m thick) superficial zone of high velocity gradient, 
overlying a zone of decreased velocity gradient extending to a depth of 
at least 150 to 200 m, the maximum depth accessible to the method. 

Based on the velocities of these two zones, except the bottom 
zone at site 14, the sites fall into two groups. At the older sites 
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(14 and 16), the shear velocities of the piecewise linear model and 
the average congressional velocities in the uppermost 10 to 13 meters 
are lower, and the uppermost layers of the layered models are thicker 
than at the younger sites (12 and 15). The velocities are lower at site 
16 than at the other sites. The lower velocities at the older sites 
may simply reflect longer exposure to meterold bombardment at these 
sites. 

The velocity gradients of the bottom zone at station 14 are 
significantly higher than those at the other three sites. The higher 
velocities under station 14 at depths between 50 and 100 meters may 
be those of an older surface layer preserved by deposition of the Fra 
Mauro formation 4 billion years ago, i.e., at the time of the Imbrium 
impact. Alternatively , the higher velocities may correspond to the 
basal portion of the Fra Mauro Formation Itself, indurated by sintering 
at the time of deposition. 

In conclusion, the observed horizontal-to-vertical spectral 
amplitude ratios were shown to be characteristic of each site and 
were fitted with the theoretically calculated spectral ratios. Owing 
to differences between the spectra of the two horizontal components of 
ground motion, the limitations caused by using only the fundamental 
mode Rayleigh waves for ellipticity modeling, and disregarding lateral 
inhomogeneities, this technique may not be used to model fine details 
of the structure. The results should be interpreted as representing 
gross velocity distributions as shown by the rather good fit of the 
piecewise linear shear velocity models to the obtained layered models. 



PART 4. ANALYSIS OF THE LONG-PERIOD LUNAR SEISMIC SIGNALS CORRECTED 
FOR INSTRUMENTAL RESPONSE AND NEAR-SURFACE EFFECTS. 

4.1. Introduction 

In the first two parts of this dissertation, I described how 
Inverse filters were designed to remove the Instrumental response and 
the near-surface effects from the long-period lunar seismograms. The 
application of these Inverse filters resulted In comparable seismic 
signals from the various seismographs. (See Section 2.4. and Figure 
33). It was expected that the seismic signals containing no dominant 
spectral peaks would yield better-defined direct shear and secondary 
wave arrivals than those derivable from the original, uncorrected 
seismograms. There are large uncertainties in the seismic event 
locations and in the velocity distributions (see Section 2.) because 
the small number of seismic stations (4), the small seismic sources, and 
the intensive scattering produced few seismic arrivals which can be 
observed on the lunar seismograms with uncertainty less than a few 
seconds; most first arrivals are emergent and most direct shear and 
secondary arrivals are buried in the scattered wavetrain. 

Identification of direct shear wave arrivals Is important in 
order to narrow the bounds on the velocity distributions for the lunar 
interior and on the source location. The recognition of secondary 
arrivals may clarify the depths of velocity discontinuities in the 
lunar interior. Velocity discontinuities were identified at depths 
of 20 to 25 km in the crust (Toksoz et al . , 1974; Goirs, 1978; Goins, 
_al., 1978), at depths of 45 to 90 km for the bottom of the crust 
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(Latham et al.. , 1972; Toksoz et al.. , 1974; Dainty et a].. , 1977; Goins, 
1978; Koyama and Nakamura, 1979), and at depths of 300 to 480 km In the 
mantle (Nakamura g£.gl. , 1974b and 1976; Voss e^gh , 1976; Latham 
£l£l_. , 1977; Goins, 1978; Goins §t ah , 1978). Well -constrained first- 
order discontinuities would narrow the bounds on the average velocities 
In the upper mantle. The determined values of the velocities In the 
upper mantle vary from 8.1 to 7.75 km/sec for the compresslonal waves 
and from 4.7 to 3.7 km/sec for the shear waves. The velocity gradients 
determined for the mantle also vary widely (e.g., Latham et al_. , 1973; 
Dainty fit al.. , 1974; Nakamura et al_. , 1974b and 1976; Goins et al_. , 1978; 
Koyama and Nakamura, 1979). 

To search for secondary arrivals, I located the seismic events 
with assumed velocity distributions, and then seismic record sections 
were made for the three source types, impacts, shallow moonquakes, and 
deep moonquakes. Seismic sections of the radial, transverse, and 
vertical components of ground motion were examined for consistent sets 
of secondary arrivals, as discussed in Section 4.3. The difficulties 
in locating the seismic events are presented in the following section. 

The most important question to be answered is whether the compresslonal 
and shear wave arrivals are satisfactory to distinguish among the 
possible velocity distributions in the lunar interior; primarily, 
whether they are able to resolve the thickness of the crust and the 
average mantle velocities. 

4.2. Direct compressional and shear wave arrivals, and locations of 
major lunar seismic events. 
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The sparse array, strong scattering, emergent beginnings, 
unrecognizable shear wave arrivals from surface sources, and occasional 
glitches Impose severe limitations on the number of locatable seismic 
events. Furthermore, rather large differences occur In the observed 
arrival times for stations 12 and 14, which are relatively close to 
each other. These differences are partially attributed to the different 
amplifications of the near-surface zones at stations 12 and 14, and to 
the facts that the short-period vertical seismometer at station 12 did 
not work normally and that the long-period vertical seismometer at 
station 14 only operated Intermittently. Thus, most of the first 
arrivals were read from the long-period vertical comoonent for 
station 12, and from either the short-period vertical component or the 
horizontal components for station 14. The number of reliable direct 
body- wave arrivals observable for a seismic event Is farther limited 
by the facts that impacts produce weak shear waves and that moonquakes 
tend to produce weak compres:ional waves. The amplitudes of the seismic 
wave arrivals, especially those from Impacts, are further obscured by 
intensive scattering. 

The emergent beginnings and obscured shear wave arrivals resulted 
in different estimates for the arrival tires, as read by different 
investigators (see Table 20). The observed-minus-calculated arrival 
times, the residuals, may not be distributed randomly for emergent 
arrivals. This contradicts the assumptions involved in finding the 
locations of the seismic events by minimizing the residuals in the 
least-squares sense (e.g., James et a]_. , 1969). To Increase the 
reliability of the observed arrival times, I read the direct shear wave 
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arrivals from the long-period seismic records and Ignored the emergent 
first arrivals when better-defined arrivals could be observed at other 
stations. This should lead to Increased stability of the least- 
squares solutions of the locations (Ibid.). The direct shear arrivals 
were read using both the original seismograms and the seismic records 
corrected for the Instrumental response and near-surface effects. 
Preliminary seismic sections containing records filtered with a three- 
pole Butterwurth filter In various frequency bands were also used to aid 
the search for consistent sets of direct shear arrivals. Generally, It 
was easietto read the direct shear arrivals from deconvolved seismic 
records filtered In a frequency band of 0.3 to 0.5 Hz and rotated Into 
the transverse direction. The seismic sections and the rotation of the 
horizontal components of ground motion are described In the following 
section. 

As described In Section 2.4., from the complete data set I 
Initially selected 65 large events (see Table 9). I expected to obtain 
the locations of these events with residuals less than a few seconds. 

Two different velocity models were used to locate these events. These 
two models are given In Tables 16 and 17 and are shown In Figures 44 and 
45, Velocities In Table 16 correspond to those published by Nakamura 
et al. , (1976); velocities in Table 17 are similar to those of Koyama 
and Nakamura (1979), Goins aK (1978), and Goins (1978). These 
velocity models were only working models; they were used to locate the 
seismic events and to calculate theoretical arrival times *or reflected 
and converted seismic waves. 


As seen from Figures 44 and 45, the velocities are continuous. 
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except at depths of 25, 60 or 45, 200, 425 and 1350 km. The continuous 
velocities were obtained by Mohorovlclc's relation (Bullen, 1965, 
Chapter 7): 


v(r) « v^(r/r^ ) a 


(25A) 


where v Is the velocity Inside a spherical shell at radius r, v^ Is the 
velocity at the outside radius of the shell, r^. If the velocity is 
Vj +1 at a radius r*r^ + ^ (where r^ar^), then 


a 


ln(v 1+1 /v ; 1 


(356) 


Tables 16 and 17 specify two congressional and two shear velocities at 
certain depths. The velocities change abruptly at such depths, which 
represent first-order discontinuities, where reflection and/or 
conversion can take place. 

If the angle of Incidence is 1. at radius r • , then the radius of 
the deepest point of the downgoing ray, Rj, is 

Rj = r.(s1n Ij) 1 ^ 1- *) (35C) 

If the seismic ray traverses this shell only once, then the contributions 
of this shell to the travel time and distance, T. and 0^, respectively. 


are 
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- 7^* cos 1 i+1 , 

v 1+l 1+1 J 




i r 4 



t 


If R i±r u] (350) 
If (35E) 

If R 1 > r i+1 (35F) 
If R 1 ± r 1+1 (35G j 


where 1^ Is calculated from the following relation 

r 1 s1n *1 Vi s1n Vi 
V 1 " Vl 


( 35H) 


The travel times and distances for the upgolng rays are given by (35E) 
and (35G), respectively, and 1^ Is calculated from (35H) since *1+1 is 
known. 

Some of the weak Impacts listed In Table 9 could not be located 
without having large residuals. Table 18 lists those events which could 
be located with residuals of a few seconds or less. The origin times 
and locations In Table 18 are based upon the observed arrival times In 
Table 19 and upon the velocity distributions In Table 17. Table 19 also 
lists the calculated arrival times. I assumed surface sources for events 
classified as Impacts and a source depth of 80 km for shallow moonquakes 



81 


because the focal depths could not be resolved (e.g., Nakamura et al.. , 
1979). Figures 46 and 47 give the epicenters of the Impacts and of 
the shallow moonquakes, respectively. I also attempted to locate the 
deep moonquakes, but the observed arrival times led to residuals of 
up to 10 to 15 seconds. Therefore, I did not list the observed arrival 
times for the congressional waves In Table 19, and I decided to use 
those deep-moonquake locations which were published by Lammleln (1977). 
Although Lammleln 1 ; locations may not be consistent with the observed 
shear wave arrival times In Table 19, the seismic record sections 
for deep moonquakes may be used to search only for those seismic wave 
arrivals which are converted at or reflected from a discontinuity 
closer to the surface than to the source. 

I also decided to use events H07, H08, H09, C02, and C03, which 
were recorded at two stations only. For these events, the locations in 
Table 18 are preferred to those in the footnote of Table 18. The 
seismograms of HOI and H07, as well as of C01 and C02, are very similar 
at a given station. The first several seconds of the seismograms of 
C03 and C04 also show some degree of similarity. Finally, the loca- 
tions for H08 and H09 are near mare areas like those for other shallow 
moonquakes located in this study. The locations and origin times for 
the artificial impacts were determined by spacecraft navigation and are 
given by Lauderdale and Eichelman (1974), except for event CS6, which 
had to be located by selsmological methods because the tracking signal 
was lost before impact. 

The origin times, locations, and residuals calculated for 
natural impacts and shallow moonquakes with the velocities in Table 16 


C-T- 
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(source depth of 100 km was assumed for shallow moonquakes) are 
very similar to those calculated with the velocities In Table 17. 

Since most residuals ahe within the uncertainties In the observations, 
the arrival times of the direct compress lonal and shear waves from 
natural lunar seismic sources cannot be used to distinguish between 
these two velocity models. The uncertainties In the observations were 
estimated by the differences In the arrival times as read by the 
different investigators (see Table 20). For artificial Impacts from 
distances between 11.2 and 35.7°, the residuals calculated with the 
velocities in Table 17 are a few seconds smaller than those calculated 
with the velocities in Table 16. Thus, the distances for aligning the 
seismic records in the sections and the arrival times for reflected and 
converted seismic waves were calculated with the velocities in Table 17, 
as discussed in the following section. 

One important result of the locations obtained in this study is 
that two shallow moonquakes, HOI and H07, occured at the same source 
region almost two years apart. This finding provides further evidence 
that shallow moonquakes are true tectonic quakes (Nakamura et al_. , 1979). 

4.3. Secondary arrivals. 

It was expected that secondary arrivals would be recognized 
most readily when the true ground motions could be compared among the 
various sites. Thus, the seismograms were corrected for instrumental 
response and near-surface effects, were filtered and were rotated into 
the radial, R, and transverse, T, directions according to 
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R(t) * -X(t)cosa - Y(t)s1n<a (36A) 

T(t) » X(t)s1na - Y(t)cos a (36B) 

where X(t) and Y(t) are the deconvolved and filtered X and Y components 
of ground motion as the function of time, t, and a Is the angle from 
the positive X axis of the seismometer to the direction toward the 
source from the station (positive Is clockwise), a's are listed In 
Table 18 under the column 'Azimuth'. R is positive away from the 
source, T is positive to the right looking from the source. 

Seismic sections were made for the three major source types, 

artificial, shallow moonquakes, and deep moonquakes, by aligning the 

radial, transverse, and vertical components of ground motion according 

to the distance between the source and the station (see Table 18), 

and using the arrival times of either the direct compress ional or the 

direct shear wave (see Table 19). The seismic record sections given in 

the pocket consist of seismic records filtered with a three-pole 

o 

Butterworth filter (Bath, 1973, p. 260) having high-pass and low-pass 
frequency cut-offs at 0.4 and 1.5' Hz, respectively. The signal ampli- 
tudes on these seismic record sections are positive toward 
increasing distances. The seismic records are identified by a 
sequence number, a four-character word, and a scale factor used for 
scaling the amplitudes of the seismic records by division (see also 
Seismic Section Captions). 

If one of the components was not working normally (see Table 9), 
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then the corresponding seismic record was zeroed. This process 
resulted In plotting one horizontal component of ground motion on 
both the radial and transverse sections when the other one was 
zeroed. No seismic record was plotted when either both horizontal 
components or the vertical component were missing. The glitches were 
also zeroed, and the seismic record was not plotted at those times when 
the amplitudes were clipped. Since some of the seismic records are 
shifted slightly along the distance axis for ciarlty of the sections, 
a line is drawn between the correct distance (see Table 18) and the 
beginning of the record, which is marked by a small tlckmark. The 
background noise may be estimated for Impact and shallow-moonqvake 
seismic records from that portion of the record which lies between 
this tickmark and the compresslonai wav arrival. 

Other seismic record sections, not included in the pocket, were 
also made using cut-off frequencies other than 0.4 and 1.5 Hz In order 
to compare the seismic signals in various frequency bands. The lengths 
of the seismic records were also varied In order to compare various 
parts cf the seismic records. Polarization filters were also used for 
enhancing rectilinear motions, as discussed below. 

Low coherency among the displacement components of ground 
motion and weak onsets of seismic arrivals In comparison with the 
amplitudes of the ambient signal make the Identification of direct 
shear and secondary arrivals difficult. To enhance the rectilinear 
motion, thus to reduce the number of possible readings for direct shear 
and secondary wave arrivals, polarization filters were designed by 
Voss et aK, (1976), Jarosch, (1977), Dainty etaj_. (1977), Goins (1978), 
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and Goins et aK (1978). The practical and theoretical difficulties 
In applying the polarization filters for the long-period lunar seismo- 
grams are significant. In addition to the low coherency among the 
various ground motion components and to the strongly scattered body 
wave arrivals, unknown relative magnification among the seismographs, 
almost vertically arriving body waves, missing or very small vertical 
components of ground motion, and the fact that the polarization filters 
also pass signals which are energetic only on one component made the 
application of the polarization filters doubtful. In fact, the polari- 
zation filters did not aid the identification of secondary arrivals, 
because most arrival times read from the polarization-filtered seismic 
record sections coincided with the increase in the signal amplitudes 
only on one component of the ground motion. 

Apart from the increase of the signal amplitudes, no criterion 
could be set to read any set of secondary arrivals. Theoretical travel 
times were calculated for many possible reflected and converted seismic 
waves and were compared with the times where the signal amplitudes 
Increase. Uncertainties in the distances and in the arrival times at 
which the seismic records are aligned were estimated to result in 
seismic waves arriving up to several seconds earlier or later than pre- 
dicted by theoretical calculations. Possible variations in the veloci- 
ties and the fact that the source depths for moonquakes cannot be deter- 
mined accurately additionally complicate the identification of the 
converted and reflected seismic wave arrivals. 

The most energetic reflections are normally those from the free 
surface, which may be important in the impact and shallow moonquake 
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seismograms. The theoretical arrival times for shear waves reflected 
once, twice, and thrice from the free surface (SS, 3S, and 4$, respec- 
tively) were calculated with the model In Table 17 and are given on the 
Impact sections. No consistent Increases In the signal amplitudes 
were found at arrival times predicted for compresslonal waves reflected 
from the free surface. The most Impulsive arrivals observed for SS are 
at distances 4 to 5.3°, for 3S are between 6.5 and 9°, and for 4S are 
between 9 and 12° (see the transverse component of ground motion). 

These arrivals Indicate that the direct shear wave arrivals are 
rather energetic between 2 and 3°, which, according to models in 
Tables 16 and 17, corresponds to an abrupt Increase In the shear velo- 
cities at a depth between 20 and 30 km (see Sullen, 1965, Chapter 7). 
Possible peg-leg shear wave multiples In the upper crust, corresponding 
to the theoretical travel time curve S25, also suggest a velocity 
discontinuity at a depth of about 25 km. Peg-leg shear wave multiples 
in the whole crust, corresponding to the theoretical travel time curve 
S45, appear to be recognizable at distances of 70 to 90°. The direct 
shear waves are also slightly more impulsive at distances greater than 
70° than between 50 and 70°. The increased Impulsiveness of the shear 
wave arrivals at distances of 70 to 90° indicates that the shear wave 
velocities decrease at a smaller rate below 200 km than above this 
depth. The absence of recognizable shear waves at distances greater 
than 90° may be the result of either Increasing attenuation or an 
Increasingly negative shear wave' velocity gradient at depths from 400 
to 450 km. 

Although the minimum distance at which surface reflections 
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can be observed Increases with Increasing source depth, the SS arrivals 
do not constrain the depths of shallow moonquakes because no consistent 
SS arrivals were recognized on the shallow-moonquake sections. The 
amplitudes of some of the shall ow-moonquake signals Increase at arrival 
times predicted for peg-leg multiples In the upper crust and In the 
whole crust (see curves S25 and S45 on the transverse component). 

S to P conversions at the 25 and 45 km discontinuities can also 
be observed at about 5 and 8 seconds before the direct shear arrival on 
some of the vertical seismic records for deep moonquakes. Latham et al . 
(1970) explained high frequency trains of waves beginning on the 
horizontal components at the Apollo 14 site about 8 sec after the direct 
compressional waves as shear waves converted at a depth of approximately 
25 km. Based upon models of Tables 16 and 17 and upon similar arrivals 
at 4.5 to 5 seconds after the first arrivals, I Interpret the arrivals 
following the first arrivals and preceding the shear wave arrivals by 
4 to 5 and 8 to 9 seconds as seismic waves converted at depths of 
approximately 25 and 45 km, respectively. The depths of discontinuities 
where these conversions take place cannot be constrained by more than 
3 and 5 km, respectively, because the uncertainties of these observations 
are about 1 second. 

Seismic arrivals preceding the direct shear wave arrivals by 
20 to 30 seconds were observed by Latham et ajL (1974) and by Nakamura 
et aU (1974b) on the vertical seismograms from deep moonquakes. Accord- 
ing to velocity distributions in Table 17, these seismic arrivals can be 
identified as conversion of the shear wave to a compressional wave at a 
depth slightly greater than 200 km (see the theoretical arrival time 
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curve for the shear-to-compresslonal wave conversion at the 200 km 
discontinuity on the deep-moonquake seismic record sections-) 

Beside the high-frequency trains following the first arrivals 
by about 8 seconds* which are most characteristic at site 14* no 
station dependency of any of the above discussed secondary arrivals was 
found. 


4.4. Surnnarv of Part 4. 

In this part of the dissertation-, the study of the long-period 
lunar seismograms corrected for Instrumental response and near-surface ef- 
fects has been described. It was my objective to Identify consistent sets 
of direct shear and secondary wave arrivals, and thus to constrain the 
velocities In the lunar mantle and the depths of the velocity discon- 
tinuities. Two velocity models, one with a thicker crust (55 km) and 
higher upper mantle velocities (8.1 and 4.6 km/sec) and the other with a 
thinner crust (45 km) and lower upper mantle velocities (7.7 and 4.4 
km/sec), were used to locate the natural Impacts and the shallow moon- 
quakes and to obtain the travel time residuals. Although most residuals 
for both models are within the uncertainty of the observed arrival 
times, the model with thinner crust and lower upper mantle velocities 
Is preferred because It gives smaller residuals for the arrival times 
observed for artificial Impacts at distances between 11° and 36°. Seismic 
sections were made of the radial, transverse, and vertical components of 
ground motion In order to search for consistent sets of secondary 
arrivals. The only criterion for reading the secondary arrivals was the 
Increase of the signal amplitudes. The velocities and the depths of the 
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velocity discontinuities cannot be severely constrained by 
secondary arrivals. However, based upon the amplitudes of the shear 
wave arrivals reflected from the free surface and upon the arrival 
times of converted waves, It Is likely that the crust Is thinner than 
55 km and that a first-order discontinuity separates the upper and lower 
crust at a depth of about 25 km. Also, some Indications were found for 
a first- ora second-order discontinuity at a depth of about 200 km. 
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CONCLUSIONS 

7h« small number of lunar seismic stations (4), the small 
magnitudes of seismic sources, and strong scattering Imposed severe 
limitations on locating the natural seismic events and on deriving 
structural information about the lunar Interior from the available seis- 
mic data. Since these conditions cannot be Improved for the time being. 
It was necessary to correct for those effects which Influenced the 
detection of seismic signals. In this study, Inverse filters were 
designed to correct for the Instrumental response, for the coupling of 
the seismometer to the ground, and for the near-surface structural 
effects. By removing the predictable digitizing errors and by applying 
these Inverse filters, the calculated ground motions were found to be 
more comparable at the various recording sites than were the recorded 
ground motions. The deconvolved signals permit reading the shear wave 
arrivals with smaller uncertainty than Is possible from the original 
seismograms, and thus yield better-constrained locations. However, no 
conclusive set of secondary arrivals could be recognized on the seismic 
sections made for Impact, shallow-moonquake, and deep-moonquake records. 
Although this study of the long-period seismic signals did not result In 
severe constraints on the structure of the lunar Interior, the transfer 
function of the long-period seismographs and the near-surface 
velocities were derived from the analysis of the calibration 
pulses and of the horlzontal-to-vertlcal spectral amplitude ratios of 
the long-period lunar seismograms, respectively. 

The finding of consistent sets of secondary arrivals remains the 
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most Important goal for the analysis of the lunar seismic signals In the 
future. Beside further limiting the number of Impacts and shallow 
moonquakes to only the largest of those used In this study t stacked 
deep-moonquake seismic signals should be used as the continuation of 
this work. 



APPtritflX: LEAST-SQUARES FIT TO MON-LINEAR FUNCTION WITH A 
LINEARIZATION OF THE FITTING FUNCTION. 


If the observations y^ at x^ (1*1 , 2,..., M) are to be fitted 
with a function y(x,aj»a 2 »»*.»a n ), which Is the function of the Inde- 
pendent variable x and of the parameters a^(j«1»2,...,n) then the 
method of linear least squares can be used to determine the optimum 
values for the parameter Increments 6a First, the fitting function 
Is expanded to first order In a Taylor's expansion as a function of 
the parameters: 

n 3>' ft (x) 

y(x) » y 0 (x) ♦ i - 5 j— 6a j (Al) 

j-1 J 

where y Q (x) - y(x,aj ajj). 


and 




are the starting values for the set of parameters. 


The partial derivatives, y Q j(x), may be calculated either analytically 
or numerically from the analytical expression of the fitting function. 
The differences between the observed and calculated values are 

given by 
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The last equation can be abbreviated by a matrix equation: 


0 - Ax ■ £ (A2A) 

where once-underlined lower-case letters Indicate column vectors, and 
twice-underlined capital letter denotes * matrix: 




Is the observed-mi nus-calculated column vector. 


| y ol (x 1> ••• y on (x l> 


y ol^ x H> •" J'oi/’W 


Is the matrix of the partial derivatives. 


5*i 


«a. 


(A2B) 


(A2C) 


(A2D) 


Is the parameter increment or adjustment column vector, and 
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3 

r l 

• 
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V y<X M> 


? M 


(A2E) 


is the residue or error column vector. 

The unknown vector, x, is calculated by minimizing the sum of 
the squared residues, that is the mean-square error, E, 


e * linn 


M 

Z 

1*1 




[b-Ax || = minimum 


(A3) 


M must be equal to or greater than n, and the column vectors of matrix 
A must be linearly independent. This process can be repeated with the 
new values of the parameters (a.+6a., j*l,2,...,n) as long as the 

J J 

Iteration converges, that is as long as the mean-square error decreases. 
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Table 1. The lunar seismic network.^ 

Station Installation Location (degrees) Azimuths^ Total number of days^ 

Latitude Longitude (degrees) In peaked In fTat mode 


11 

July 

21,1969 

0.68 

N 

23.42 

E 

0 . 


21' 

12 

Nov. 

19,1969 

3.04 

S 

23.42 

W 

180. 

2055 7) 

81 T 

14 

Feb. 

5,1971 

3.65 

S 

17.48 

W 

0 . 

2288 8 * 


15 

July 

31,1971 

26.08 

N 

3.66 

E 

0 . 

1615 

638 1 

16 

Apr. 

21,1972 

8.97 

S 

15.51 

E 

-25.5 

1350 

638* 


1) The network was turned off on September 30, 1977. 

2) The listed azimuths are the directions of the ground motion which give 
positive signals on one of the horizontal components (X) and no signal 
on the other (Y). Azimuths are measured clockwise from North. The 
direction of the ground motion which gives positive signals on Y Is 
90° clockwise from the direction listed. 

3) Only those days were excluded from this count when the signal was lost 
for at least one full day. This happened at station 14 only between 
March 1 and March 5, 1975, between January 18 and February 19, 1976, 
between March 17 and May 20, 1976, between June 8 and June 10, 1976, 
between October 9 and November 12, 1976, and between July 30 and 
August 4, 1977. 

4) Station 11 was turned off for the lunar night between August 2 and 
August 17, 1969. 

5) Station 12 was operating in the flat mode between November 19 and 
November 22, 1969, and between October 14, 1974 and April 9, 1975. 


preceding page blank NOT 
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Table 1. (continued) 

6) Stations 12, 15, and 16 were operating In the flat mode between 
June 28, 1975 and March 27, 1977. 

7) The short-period vertical seismometer at station 12 never worked 
properly. 

8) The long-period vertical seismometer at station 14 did not work 
properly most of the time. 
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Table 2. Nominal values of the seismograph constants. 


Constant 

Nominal value 

Parameter 

A/0 converter transfer function 

K • 0.2049 DC/-*' 


Displacement transducer transfer f. 

K<| ■ 5000V/cm 


Output gain 

K 3 ■ 31.6 at 0 dB 1 2) 

*1 " ¥3 

Pendulum damping constant 

0 » 0.85 

a 2 « 2s 

Pendulum free period 

2ir/u» * 15 sec 

a, » u 


0 

3 0 

Demodulator low pass cutoff frequ. 

u d » 47.62 sec" 1 

a 4 ’“d 

Output high pass cutoff period 

2w/u»_ * 100 sec 

a 

a 5 " “a 

Output low pass cutoff period 

2ir/wj ■ 0.72 sec 

a g • oi 1 

Coll -magnet transfer function 

K 2 - 0.0016 gal/V 

a 7 “ K 1 K 2 

Feedback low pass cutoff period 

2 ir/u f * 6300 sec 

a 8*“f 


1) As referred to In the text and In the figures. 

2) The attenuator had four settings with lOdB increments from 0 dB 
to 30 dB (0 dB was used most of the time) . 
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Table 3. 


Description of the calibration pulses. 


1 ) 



Peaked mode^ 




Flat mode 




(‘74, 

day 175) 




(1976. 

day 149) 



Comp 

Maxima Length^ 

Maxima 

Minima 

Zero 

length^ 

Noise 








crossing 

(DU's) 


(DU's) (sec) 

(DU's) (sec) 

(DU's) 

(sec) 

(sec) 

(sec) 

12X 

14.5 

45 

257.0 

7.6 

-49.0 

43.0 

24.0 

150 

15 

12Y 

14.5 

42 

236.0 

7.0 

-40.5 

42.0 

22.5 

60 

28 

122 

15.5 

43 

361.0 

9.8 

-189.5 

34.3 

22.2 

11C 

9 

14X 

14.5 

35 



Station 

14 Is 

unstable 



14Y 

16.0 

35 








142 4 > 





In flat mode 



14.5 

45 








15X 

15.5 

40 

281.5 

7.6 

-54.5 

45.0 

23.1 

160 

7 

15Y 

15.5 

40 

304.0 

8.0 

-67.0 

39.4 

21.5 

120 

11 

15Z 

124.0 

80 1 2 * * 5) 

511. 5 6) 

-345.0 

44.0 

23.6 

210 

6 

16X 

19.5 

45 

318.0 

7.4 

-64.0 

41.0 

23.4 

no 

11 

16Y 

13.5 

40 

263.5 

7.7 

-56.0 

40.0 

21.2 

60 

16 

16Z 

19.0 

50 

338.5 

7.5 

-71.0 

38.0 

20.6 

100 

11 


1) Quantities were measured from the average of the positive and the sign- 
reversed negative calibration pulses (see Figures 3 and 4). 

2) The noise was about 1 digital unit (DU) for all calibration pulses In 

the peaked mode, except for component 16Y and for the positive calibra- 

tion Dulse for component 15Z, for which the noise '"as 2 DU's. 

^ Length refers to the approximate duration of the calibration pulse. 
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Table 3. (continued) 

that Is the time difference between the beginning of the calibra- 
tion pulse and the return to the background noise level. 

41 

' Calibration pulse of 1971 day 170 was used for component 14Z. 

The positive calibration pulse follows 56 seconds after the 
negative one, while the signal from the negative pulse Is still 
about 2 DU's above the noise level. 

^ The amplitudes of the calibration pulse are clipped. 
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Table 4. 


Seismograph constants obtained in this study. 


1 ) 



Peaked mode (2ir/»«.*15sec) 

Op 

Flat mode (2w/#y«6300sec) 

# f ■ 

Comp 

C P 

s p 

*7p 

Peak 2 * 

c f B f le. *7f 

“bf 

!ze 

a 7f 

“of °P 




gal /cm freq. period 

1 sec gal/cm 


12X 

0.92 

0.726 

7.04 

0.430 

2.32 

0.85 0.670 14.9 6.06 

1.16 

0.731 

12Y 

0.95 

0.763 

7.50 

0.443 

2.26 

0.81 0.602 14.2 6.43 

1.17 

0.723 

12Z 

0.93 

0.704 

6.87 

0.424 

2.36 

0.87 1.17 21.6 7.43 

0.92 

1.01 

14X 

0.97 

0.742 

7 73 

0.450 

2.22 

Station 14 1$ unstable 



14Y 

1.10 

0.739 

8.48 

0.469 

2.13 

In flat mode 



742 

0.93 

0.719 

7.69 

0.450 

2.22 




15X 

1.0 

0.715 

7.61 

0.446 

2.24 

1.0 0.691 14.5 6.49 

1.17 

0.692 

15Y 

1.03 

0.678 

8.10 

0.459 

2.18 

0.99 0.732 15.6 7.31 

1.11 

0.705 

15Z 

7.89 

0.670 

7.68 

0.446 

2.24 

cal pulse Is clipped 



16X 

1.13 

0.717 

6.46 

0.414 

2.42 

1.04 0.721 15*. 4 5.94 

1.09 

0.736 

16Y 

0.91 

0.726 

8.12 

0.459 

2.18 

0.85 0.691 15.4 7.36 

1.10 

0.742 

15Z 

1.22 

0.669 

7.46 

0.440 

2.27 

1.120.68015.2 7.84 

0.95 

0.678 

Median 

0.717- 

- 7.61- 



0.691 4) 15.2 4) 6.49 4 * 

1.10- 

0.723- 



0.719 

7.68 

0.446 

2.24 

1.11 

0.731 

Average 





0.94 0.684 4) 15.0 4) 6.78 4) 




1.01 3) 0.714 

7.56 

0.444 

2.25 

1.08 

0.752 


1) Table 2 lists the symbols used; C's are relative amplifications, 
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Table 4. (continued) 

C * (ca^/Cca^ for component 15X) Csee equation 00)). Subscripts 
p and f refer to the peaked mode and the fiat mode, respectively. 
Nominal values for u a » w d , and u-| are given In Table 2. 

Cstermiend from the amplitude responses to an Impulse of displace- 
ment; frequency is in Hz and period Is in seconds. 

31 

' Constant of 15Z was not used for calculating this value. 

^ Constant of 12Z was not used for calculating this value. 



Table 5. The errors in the determination of the seismograph constants 


no 

1) 



Peaked mode (X) 


Flat mode (X) 


W 

1-K£ 

C -C- 

D f 


S 

6 p 

*7p 

c f 

6 f 

at 

V 

a 7p- 

9 f' 

*7 

L f 

12X 

1 

3 

1 

6 

3 

3 

6 

-0.08 

-0.16 

0.08 

12Y 

1 

1 

1 

3 

1 

2 

5 

-0.17 

-0.17 

0.17 

12Z 

1 

5 

1 

1 

1 

1 

3 

0.06 

0.08 

0.07 

14X 

1 

3 

2 








14Y 

6 

3 

1 








14Z 

2 

6 

1 








15X 

5 

4 

1 

1 

1 

1 

2 

-0.01 

-0.17 

0 

15Y 

3 

1 

1 

1 

1 

2 

3 

-0.04 

-0.11 

0.04 

15Z 

1 

1 

1 








16X 

1 

2 

1 

7 

4 

3 

10 

-0.02 

-0.09 

0.09 

16Y 

3 

3 

2 

2 

4 

1 

8 

-0.07 

-0.10 

0.07 

16Z 

3 

6 

1 

1 

1 

1 

2 

0.01 

0.05 

0.09 


^Tables 2 and 4 give the symbols used. 





Table 6. Absolute values of the relative differences between the amplitude responses of Figure 12 and 
component 15Z in the peaked mode. 


Scale 1 ^ Frequency (Hz) 


Comp 

factor 

0.2 

0.3 

0.4 

0.45 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

1.1 

1.2 

12X 

0.33 

33.37 

38.66 

51.69 

7.40 

7.51 

8.70 

13.45 

15.75 

16.94 

17.73 

18.22 

18.58 

12Y 

0.82 

25.23 

26.02 

23.99 

0.00 

7.00 

17.25 

19.38 

20.29 

20.74 

21.03 

21,21 

21.33 

12Z 

0.85 

34.43 

41.78 

62.88 

12.51 

13.44 

4.08 

9.53 

12.22 

13.63 

14.58 

15.17 

15.60 

14X 

0.87 

13.69 

13.00 

8.05 

0.04 

10.83 

14.54 

14.74 

14,76 

14,76 

14.74 

14,74 

14,73 

14Y 

0.95 

6.52 

10.51 

24.16 

18.56 

35.32 

18.78 

12.73 

10.03 

8,68 

7.81 

7,27 

6,89 

14Z 

0.91 

9.92 

9.62 

6.79 

0.00 

6.73 

9.67 

9.97 

10.06 

10,10 

10.12 

10.13 

10.14 

15X 

0.91 

11.29 

11.53 

10.54 

0.00 

3.87 

8.37 

9.21 

9.55 

9.72 

9.83 

9.89 

9.94 

15Y 

1.03 

9.06 

11.17 

18.46 

7.71 

14.59 

3.52 

0,59 

0.73 

1.39 

1.82 

2.06 

2.27 

16X 

0.79 

55.56 

69.67 

98.12 

21.68 

18.19 

5.62 

13.65 

17,75 

19.94 

21.42 

22.36 

23.05 

16Y 

0.95 

1.03 

3.58 

13.14 

7.48 

20.62 

12.72 

9.77 

8.38 

7.68 

7.22 

6.93 

6.73 

16Z 

0.98 

6.03 

7.48 

12.79 

1.03 

5,93 

0.97 

0.45 

1.12 

1.46 

1.69 

1.84 

1.94 


^The scale factors are the ratios of the maximum of the amplitude responses of each component (see 


Figure 12) to that of component 15Z. The amplitude responses of Figure 12 were divided by the scale 


factors . 





Table 7. Absolute values of the relative differences between the amplitude responses of Figure 12 and 
component 15X in the peaked mode. 


Scaled Frequency (Hz) 


Comp 

factor 

0.2 

0.3 

0.4 

0.45 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

1.1 

1.2 

12X 

0.92 

19.54 

24.44 

36.54 

7.70 

10.15 

0.02 

3.62 

5.31 

6.25 

6.87 

7.24 

7.53 

12Y 

0.91 

12.24 

12.74 

11.57 

0.26 

3.42 

7.91 

9.05 

9,52 

9.76 

9.91 

10.01 

10.08 

12Z 

0.94 

20.74 

27.68 

46.78 

. 12.64 

15.59 

4.11 

0.24 

2.28 

3.44 

4.17 

4.66 

5.02 

14X 

0.96 

1.78 

0.86 

2.82 

0.40 

6.39 

5.31 

4.68 

4.38 

4.21 

4.10 

4.03 

3.98 

14Y 

1.05 

16.23 

20.31 

32.03 

18.70 

29.38 

9.41 

2.93 

0.26 

1.15 

2.01 

2.57 

2.97 

14Z 

1.00 

1.32 

1.84 

3.54 

0.08 

2.54 

1.11 

0.61 

0.38 

0.26 

0.18 

0.13 

0.09 

15Y 

1.14 

18.33 

20.57 

26.41 

7.71 

8.70 

4.48 

7,96 

9.39 

10.14 

10.60 

10.90 

11.12 

15Z 

1.10 

10.15 

10.35 

9.34 

0.00 

4.34 

7.73 

8.44 

8.72 

8.86 

8.95 

9.00 

9.04 

16X 

0.87 

40.02 

53.61 

75.64 

21.68 

19.48 

2.54 

4.21 

7.48 

9.36 

10.56 

11.37 

11.95 

16Y 

1.04 

11.12 

13.81 

21.75 

7.48 

15.24 

4.01 

0.44 

1.07 

1.88 

2.37 

2.70 

2.93 

16Z 

1.08 

4.70 

3.50 

2.44 

1.03 

9.53 

8.62 

8.01 

7.70 

7,52 

7.41 

7.33 

7.27 


^The scale factors are the ratios of the maximum of the amplitude responses of each component (see 
Figure 12) to that of component 15X. The amplitude responses of Figure 12 were divided by the scale 


factors 
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Table 8. Relative noise for Instrumental deconvolution In the peaked 
mode. 


Component 

Freq. (Hz) where the s/n Is 

Relative nolse 1 ^ 


2.5 1.0 

0.5 

Maximum Average 

all but 15Z 

1. 1-1.3 1.4-1 .5 

1.6-1. 7 

0.10-0.18 0.02-0.03 

15Z 

1.5 1.7 

1.9 

0.04 0.01 


1) The maximum or the average of the residuals between the amplitude 
spectrum of the calculated and observed calibration pulses In the 


frequency band of 0 to co^ divided by the amplitude spectrum of the 
calculated calibration pulse at the peak frequency (see Table 4). 
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Table 9. 


List of the most energetic long-period lunar seismic events 


1) 


Event Station 12 Station 14 Station 15 Station 16 
Yr:day Code Com WL Hr:mn:sec Com WL Hr:mn:sec Com ML Hr:mn:sec Com ML Hr:mn:sc 


71:051 A01 

XYZ 10 

15:07 

XYZ 10 15:07 

no records 

no records 

75:086 AO 1 

flat 

18:47 

1 3 



75:113 A01 

Z 3 

11:54 

1 3 11:54 

XYZ 3 11:55 

XYZ 3 11:55 

75:304 A01 

flat 

6:34 

1 3 6:34 

flat . 6:37 

flat 6:35 

73:148 A14 

weak 


weak 

weak 

XYZ 10 18:56 

72:145 A16 

weak 


weak 

weak 

Z 10 9:24 

73:001 A17 

XYZ 10 

11:34 

XY3 10 11:34 

XYZ 10 11:32 

XYZ 3 11:34 

74:116 A20 

XYZ 10 

9:20 

XY3 10 9:20 

XYZ 10 9:19 

XYZ 3 9:19 

72:341 A33 

weak 


weak 

weak 

XYZ 3 13:29 






XYZ 3 13:32 

72:138 A34 

weak 


weak 

XYZ 10 23:33 

XYZ 10 23:33 

74:166 A42 

1 Z 10 

13:21 

XY3 10 13:21 

XYZ 10 13:22 

XYZ 10 13:23 

74:084 A44 

XYZ 10 

14:33 

XY3 10 14:32 

XYZ 10 14:31 

XYZ 10 14:32 

73:243 A46 

XYZ 10 

0:04 

XY3 10 0:04 



74:343 A46 




XYZ 10 23:49 

XYZ 10 23:46 

77:107 A73 

weak 


weak 

weak 

XYZ 10 17:13 

69:324 CL2 

flat 

22:17 

no records 

no records 

no records 

71 :038 CL4 

Z 3 

0:45 

10 0:45 

no records 

no records 

71:215 CL5 

weak 


weak 

XYZ 3 3:03 

no records 

70:105 CS3 

3 

1:09 

no records 

no records 

no records 

71 :035 CS4 

XYZ 3 

7:41 

no records 

no records 

no records 

71:210 CS5 

XYZ 3 

20:59 

Z 3 20:59 

no records 

no records 
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Table 9. (continued) 


Event Station 12 Station 14 Station 15 Station 16 


Yr:day 

Code 

Com 

WL 

Hr:mn:sec 

Com 

ML 

Hr:mn:sec 

Com 

ML Hr:mn:sec 

Com ML 

Hr:mn:sc 

72:110 

CS6 


3 

21:02 

XY3 

3 

21:02 

XYZ 

3 21 :04 

no records 

72:345 

CS7 

XYZ 

3 

20:33 

3 

3 

20:33 

XYZ 

10 20:35 

X2Z 10 

20:34 

72:134 

C01 


3 

8:47 

3 

3 

8:47 

XYZ 

3 8:48 

XYZ 3 

8:48 

71:143 

C02 

XYZ 

3 

22:20 

XYZ 

3 

22:20 

no records 

no records 

71:163 

C03 

XYZ 

3 

10:51 

XYZ 

3 

10:51 

no records 

no records 

71:193 

C04 

XYZ 

3 

18:08 

XYZ 

3 

18:08 

XYZ 

3 18:08 

no records 

72:199 

C05 

XYZ 

3 

21:57:45 

XY3 

3 

21:57:45 

XYZ 

3 21 : 56:45 

Z 3 

21:57:15 



XYZ 

3 

22:01 

XY3 

3 

22:01 

XYZ 

3 22:00 

YZ 3 

22:00 

72:213 

C06 

XYZ 

10 

18:11 

3 

10 

18:10 

XYZ 

10 18:08 

3 

18:11 

72:324 

C07 

XYZ 

10 

18:26 

XY3 

10 

18:25 

XYZ 

10 18:24 

10 

18:25 

73:113 

C08 


10 

13:56 

X 3 

10 

13:56 

Z 

10 13:55 

Z 10 

13:55 

73:262 

C09 


10 

9:32 

3 

10 

9:32 

Z 

10 9:34 

XYZ 10 

9:32 











XYZ 10 

9:35 

73:269 

CIO 

XYZ 

10 

20:53 

XY3 

10 

20:50 

XYZ 

10 20:48 

X2Z 10 

20:49 

74:038 

Cll 

Z 

10 

6:21 

XY3 

10 

6:21 

XYZ 

10 6:22 

10 

6:21 

74:181 

Cl 2 

XYZ 

10 

17:44 

3 

10 

17:44 

weak 

weak 


74:187 

Cl 3 

Z 

10 

14:17 

XY3 

10 

14:17 

XYZ 

10 14:14 

Z 10 

14:16 

74:305 

Cl 4 

flat, weak 

XY3 

10 

11:44 

Z 

10 11:45 

Z 10 

11:43 











Z 10 

11:46 

74:325 

C16 

flat 

13.18 

XY3 

10 

13:18 

XYZ 

10 13:T8 

XYZ 10 

13:15 

74:343 

Cl 7 

flat, weak 

XY3 

10 

9:32 

XYZ 

10 9:33 

XYZ 10 

9:29 

74:349 

Cl 8 

flat 

9:08 

XY3 

3 

9:07 

XYZ 

10 9:10 

XYZ 10 

9:09 
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Table 9. (continued) 


Event Station 12 Station 14 Station 15 Station 16 


Yrrday 

Code 

Com ML Hr:mn:sec 

: Com WL 

Hr:mn:sec 

Com UL 

Hr:mn:sec 

: Com WL 

Hr:mn:sc 

75:064 

C20 

flat. 

weak 

1 3 10 

21:53 

XYZ 10 

21:54 

Z 10 

21:51 

75:085 

C22 

flat. 

weak 

1 3 10 

12:04 

XYZ 10 

12:07 

Z 10 

12:06 

75:102 

C23 

XYZ 10 18:16 

1Y3 10 

18:16 

XYZ 10 

18:15 

XYZ 3 

18:14 









XYZ 3 

18:18 

75:124 

C25 

Z 

3 10:05 

1 3 3 

10:05 

XYZ 3 

10:06 

Z 3 

10:05 

77:107 

C26 

XY3 10 23:35 

XYZ 10 

23:35 

XYZ 10 

23:37 

Z 10 

23:36 

77:179 

C27 

Z 10 22:25:30 

3 10 

22:26 

XYZ 10 

22:27 

X Z 10 

22:27 

76:025 

C28 

flat- 

16:09 

123 


flat 

16:11 

flat 

16:11:45 

76:109 

C29 

flat 

8:46 

123 


flat 

8:47 

flat 

8:48 

76:121 

C30 

flat 

11:54 

123 


flat 

11:55 

flat 

11:52 

76:137 

C31 

flat. 

weak 

123 


flat, weak 

flat 

12:47 

76:240 

C32 

flat. 

weak 

weak 


flat, weak 

flat 

4:31 

77:007 

C33 

flat. 

weak 

weak 


flat, weak 

flat 

3:49 

77:153 

C34 

weak 


weak 


weak 


XYZ 10 

14:26 

76:319 

C35 

flat 

23:16 

123 


flat 

23:17 

flat 

23:18:30 

77:256 

C36 


20:33 

3 

20:34 

weak 


weak 


73:171 

H01 

XYZ 10 20:24 

XY3 10 

20:24:48 

Z 10 

20:26 

Z 10 

20:27:30 

74:192 

H02 

Z 10 0:52 

3 10 

0:52 

Z 10 

0:50 

XYZ 10 

0:50 



Z 10 0:56 

XY3 10 

0:56 

Z 10 

0:54 

XYZ 10 

0:54 

75:003 

H03 

flat 

1:46 

XY3 3 

1:46:30 

Z 3 

1:46:30 

3 

1:47 





XY3 3 

1:50:30 

Z 3 

1:50:30 

3 

1:50 

75:012 

H04 

flat 

3:19 

Y3 10 

3:19 

XYZ 10 

3:17 

XYZ 10 

3:20 
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Table 9. (continued) 


Event Station 12 Station 14 Station IS Station 16 

Yr:day Code Com WL Hr:mn:sec Com WL Hr:mn:sec Com WL Hr:mn:sec Com Wl Hr:mn:sc 


73:072 H06 

XYZ 10 

8:01 

X23 10 8:01 

12Z 10 8:02:30 

XYZ 10 8:01 


XYZ 10 

8:05 

23 10 8:05 

12Z 10 8:05 

XYZ 10 8:05 

71:192 H07 

XYZ 10 

13:27 

Z 10 13:28 

no records 

no records 

71:140 H08 

XYZ 10 

17:28 

XYZ 10 17:28 

no records 

no records 


Z TO 17:32 

XYZ 10 17:32 

no records 

no records 

71:107 H09 

XYZ 10 

7:04:45 

XYZ 10 7:04:45 

no records 

no records 


XYZ 10 

7:08:45 

XYZ 10 7:08:45 

no records 

no records 

76:066 H10 

flat 

10:15 

XY3 10 10:15 

flat 10:14 

flat 10:16 

76:068 HI! 

flat 

14:43 

XY3 10 14:43 

flat 14:45 

flat 14:43:45 

TJ 

Yr:dy;Hr:mn:sec is the beginning of the 309 minute time sections. 


Code is the event designation as used in the text, in the seismic sections, 
and in Tables 18, 19, and 20. The first letter of the code is the type of 
source; A is deep moonquake, C Is impact, H is shallow moonquake. The 
two-digit numbers are arbitrary sequence numbers for natural impacts and 
shallow moonquakes and are source locations for deep moonquakes. The 
artificial impacts are denoted by the letter C followed by a letter, 
giving the source (L is for Lunar Module and S is for Saturn IV booster), 
and by a number, giving the second digit of the mission number. WL is 
the waterlevel parameter (in *) used for Instrumental deconvolution. Com 
is 1, 2, or 3 if the X, Y, or Z component, respectively, did not operate 
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Table 9. (continued) 

properly at the time of the event, and X, Y, or Z are used to denote 
those time sections which were recorded In the peaked mode and were 
used for computation of the whitening filters, 'flat' Indicates 
that the long-period seismographs were operated In the flat mode 
at the time of the event and the seismograms of all three components 
were used In the seismic sections, 'weak* Indicates that the signal 
amplitudes were too small to be used In this study. 



119 


Table 10. Number of time section:- classified according to source type 
and to components. 


Deep moonquakes Shallow moonquakes Impacts Total 

Peaked Flat Peaked Flat Peaked Flat Peaked Flat 



Good 

Noisy 


Good 

Nol sy 


Good 

Noisy 


Good 

Noisy 


12X 

5 

2 

2 

7 

3 

4 

14 

10 

7 

26 

15 

13 

12Y 

5 

2 

2 

7 

3 

4 

14 

10 

7 

26 

15 

13 

12Z 

7 


2 

10 


4 

18 

5 

7 

35 

5 

13 

14X 

6 

1 


n 

4 


16 

7 


33 

12 


14V 

6 

3 


n 

2 


16 

11 


33 

16 


14Z 

1 



5 



5 



11 



15X 

7 


1 

1 

5 

2 

21 

3 

4 

29 

8 

7 

15Y 

7 


1 

1 

5 

2 

21 

3 

4 

29 

8 

7 

15Z 

7 


1 

8 


2 

24 


4 

39 


7 

16X 

11 

1 

1 

5 

3 

2 

11 

14 

7 

27 

18 

10 

16Y 

11 

1 

1 

5 

3 

2 

11 

12 

7 

27 

16 

10 

16Z 

12 



6 

2 

2 

22 

3 

7 

40 

5 

10 

^The number of 

time sections listed 

In 

column 

'Good" 

were 

recorded in 

the 


peaked mode and were used for computation of the whitening filters, in 
column 'Noisy' were not so used. The number of time sections listed in 
column 'Flat' were recorded in the flat mode and were used in the seismic 
record sections. 
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Table 11. Relative amplification of the horizontal components (Y/X) 1 


Sta 

Oeep moonquakes 

Shallow quakes 

Impacts 

Total 


0.2-1. 5 

0-2.5 

0.2-1 .5 

0-2.5 

0.2-1 .5 

0-2.5 

0.2-1 .5 

0-2.5 

12 

1.02 

0.99 

1.02 

0.93 

C.98 

0.91 

0.98 

0.92 

14 

1.49 

1.29 

1.14 

1.21 

1.23 

1.21 

1.21 

1.21 

15 

1.30 

1.27 

1.26 

1.13 

1.40 

1.27 

1.40 

1.26 

16 

1.39 

1.21 

1.36 

1.27 

2.07 

2.00 

1.80 

1.70 


^Average amplifications In two frequency bands, 0.2 to 1.5 Hz and 
0 to 2.5 Hz. 
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Table 12. 


Structural parameters of the near-surface zone; model 



p 

a/ 8 

Station 12 
z a 

Station 14 

Z 8 

Station 15 
z a 

Station 16 

Z 8 

1.5 

2.9 

1.2 

36 

2.9 

36 

1.7 

36 

4.4 

36 

1.6 

2.7 

4.1 

62 

5.0 

49 

6.2 

62 

9.9 

56 

1.7 

2.4 

11.3 

96 

8.6 

68 

13.0 

96 

25.1 

96 

1.8 

2.2 

19.5 

136 

18.0 

98 

25.4 

136 

49.6 

126 

1.9 

2.0 

41.1 

204 

28.1 

117 

59.3 

204 

104.6 

204 

1.9 

2.0 

84.3 

226 

64.4 

183 

104.6 

283 

168.4 

283 

1.9 

2.0 

127.5 

340 

86.4 

277 

149.9 

340 



1.95 

2.0 

•1 

430 

m 

360 

m 

480 

m 

400 


i 1 

p Is density in g/on , z Is depth In meters, a and a are congressional 


and shear velocities In m/sec 
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Table 13. 


Structural parameters of the near-surface zone; model 





Station 12 

Station 14 

Station 15 

Station 16 

p 

o/8 

z 

6 

z 

8 

z 

0 

z 

8 

1.8 

2.6 

1.9 

40 

3.2 

40 

2.8 

40 

5.1 

40 

1.85 

2.5 

7.5 

68 

5.5 

55 

9.5 

68 

10.9 

62 

1.9 

2.25 

10.3 

93 

9.5 

76 

12.8 

93 

27.1 

106 

1.9 

2.16 

18.7 

135 

19.0 

109 

22.8 

135 

51.5 

150 

1.95 

2.08 

33.6 

167 

30.3 

131 

40.7 

167 

111.2 

227 

1.95 

2.06 

80.1 

205 

70.8 

204 

74.2 

205 

178.5 

313 

2.0 

2.0 

135.9 

302 

95.5 

309 

141.2 

302 

m 

440 

2.0 

2.0 

•e 

395 

m 

400 

ee 

395 




^0 Is density In g/cm\ z Is depth In meters, a and s are congressional 
and shear velocities In m/sec. 
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Table 14. Description of the velocity distributions of model *B ' In 
the forms of 6 * a + bz and a = cz*\ ^ 


Station 

Top zone 


Bottom zone 

Power law 

N 

a 

a 

b 

h 

a 

b 

c 

P 

12 

33.7 

6.9 

16.1 

118 

1.7 

106 

0.31 

3 

159 

14 

29.6 

6.1 

10.1 

62 

2.9 

89 

0.32 

3 

134 

15 

31.7 

5.7 

19.2 

107 

1.8 

93 

0.34 

3 

156 

16 

29.8 

4.1 

25.2 

95 

1.5 

75 

0.85 

2 

126 


^ 6 and a are the shear and congressional velocities, respectively, in 
m/sec, z is depth in m, a, b, c, and p are constants, h is the depth in m 
where the gradient (b) changes its value, N is the number of uppermost 
layers in which the power law adequately describes the compressional 
velocities, and a is the average velocity in the zone consisting of these 
1 ayers . 
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Table 15. Travel times from the Apollo 14 and 16 Active Seismic 
Experiments. ^ 


X 

*m 

*1/6 

*14 

*16 

*1/3 

4.57 

55 

51.7 

44.0 

40.1 

58.3 

9.14 

91 

92.1 

87.9 

80.2 

92.6 

13.71 

123 

129.1 

131.9 

120.3 

121.5 

18.29 

151 

164.2 

175.8 

160.4 

147.1 

27.43 

206 

230.1 

245. Q 

240.6 

192.9 


^ x is the separation of the geophones in m (Gangi, 1978, Table III), t m 
is measured from the Apollo 14 and 16 stacked data in msec (ibid.), t-^g 


5/6 


is calculated for the self-compacting-powder model: t-^g 3 14.57 x 

msec (ibid.), is calculated for the Apollo 14 layered model in msec 

(ibid.), t^g is calculated for the Apollo 16 layered model in msec (ibid.), 

212 

ty 3 is calculated for the 1/3 power law model: ty 3 = 21.2 x 


msec. 
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Table 16. Velocity models used in this study with thicker lunar crust 
and higher lunar mantle velocities.^ 


Depth 

(km) 

Compr. vel. 
(km/sec) 

Shear vel . 
(km/sec) 

0.0 

0.1 

0.05 

0.25 

1.0 

0.5 

1.0 

2.1 

1.2 

5.0 

4.1 

2.4 

9.0 

5.0 

2.9 

17.0 

5.9 

3.4 

25.0 

6.4 

3.7 


6.6 

3.8 

55.0 

6.8 

3.9 


8.1 

4.6 

200.0 

8.0 

4.5 


7.9 

4.4 

425.0 

7.8 

4.2 


7.8 

3.9 

1350.0 

7.9 

3.7 


8.0 

2.0 

1 730.0 

8.0 

2.0 


^First-order discontinuities are at those depths where two compressional 
and two shear velocities are given. The velocities between the depths 
listed in this table are continuous and are given by equations (35),. 



Table 17. Velocity models used In this study with thinner lunar crust 
and lower mantle velocities.^ 


Depth 

(km) 

Compr. vel . 
(km/sec) 

Shear vel. 
(km/sec) 

0.0 

0.1 

0.05 

0.25 

1.0 

0.5 

1.0 

2.1 

1.2 

5.0 

4.1 

2.4 

9.0 

5.0 

2.9 

17.0 

5.9 

3.4 

25.0 

6.4 

3.6 


6,6 

3.9 

45.0 

6.8 

4.0 


7.7 

4.4 

200.0 

7.6 

4.35 


7.6 

4,25 

425.0 

7.6 

4.20 


7.6 

4.10 

1350.0 

7.6 

3.80 


7.6 

3.50 

1738.0 

7.6 

3.50 


^First-order discontinuities are at those depths where two compressional 
and two shear velocities are given. The velocities between the depths 
listed in this table are continuous and are given by equations (35). 



Table 18. Reference times, origin times, and locations for selected major lunar seismic events.^ 
Code Reference time Origin Epicenter (degree)*^ Distances (degrees) 3 ) Azimuths (degrees?**) 



yr:day 

hr:mn 

(sec) 

Latitude Longitude 

12 

14 

15 

16 

12 

14 15 16 

CL2 

69 

324 

22 

17 

17.7 

3.94 

S 

21.21 W 

2.4 




67.8 


CL4 

71 

38 


45 

25.7 

3.42 

S 

19.67 W 

3.7 

2.2 



-84.1 

-84.1 

CL5 

71 

215 

3 

04 

-23.0 

26.36 

N 

0.25 E 



3.1 



-84.0 

CS3 

70 

105 

1 

10 

-19.0 

2.75 

S 

27.86 U 

4.5 




93.6 


CS4 

71 

35 

7 

41 

-4.6 

8.09 

S 

26.00 W 

5.7 




26.8 


CS5 

71 

210 

20 

59 

-17.1 

1.51 

S 

11.80 W 

11.7 

6.1 



-97.3 

69.5 

CS6 

72 

110 

21 

02 

-2.4 

2.33 

N 

24.02 W 

5.4 

8.9 

35.7 


173.6 

-47.7 127.2 

CS7 

72 

345 

20 

33 

-17.7 

4.21 

S 

12.30 W 

11.2 

5.2 

34.0 

28.0 

-83.6 

96.4 -151. -56.5 

C01 

72 

134 

8 

47 

-20.0 

1.40 

N 

16.76 U 

8.0 

5.1 

31.5 

33.8 

-124. 

8.2 -138. -48.2 

C02 

71 

143 

22 

20 

-22.5 

0.65 

N 

17.23 W 

7.2 

4.3 



-121. 

3.3 

C03 

71 

163 

10 

51 

-98.5 

32.42 

N 

29.04 W 

35.9 

37.7 



171.9 

-16.1 

C04 

71 

293 

18 

08 

-108. 

30.70 

H 

28.58 U 

34.1 

36.0 

28.6 


172.1 

-16.4 -73.2 

COS 

72 

199 

21 

56 

-335. 

31.50 

N 

148.90 E 

150.6 149.3 

113.5 

131.3 

-167. 

23.2 32.0 81.0 

C06 

72 

213 

18 

08 

-4.3 

33.90 

N 

3.92 W 

41.3 

39.7 

10.2 

46.7 

-155. 

17.8 -38.2 -126. 

C07 

72 

324 

18 

24 

-160. 

61.18 

N 

16.27 W 

64.3 

65.2 

40.5 

77.4 

178.5 

-4.7 -21.7 6.3 


Table 18. (continued) 


Code Reference time Origin Epicenter (degree) Distances (degrees) Azimuths (degrees) 



yr: 

:day 

hr; 

:mn 

(sec) 

Latitude 

Longitude 

12 

14 

15 

16 

12 

14 

15 

16 

C16 

74 

325 

13 

16 

-29.3 

6.14 

S 

21.77 

E 

45.1 

39.2 

36.7 

6.8 

-84.4 

95.2 

148.9 

91.4 

C23 

75 

102 

18 

14 

-55.0 

3.20 

N 

34.37 

E 

58.2 

52.3 

37.3 

22.4 

-95.7 

83.1 

122.7 

83.4 

C25 

75 

124 

10 

05 

-362. 

35.00 

S 

131.70 

W 

103.1 

107.4 

140.9 

126.2 

52.9 

-129. 

-114. 

-121. 

C26 

77 

107 

23 

35 

-166. 

20.09 

S 

61.00 

W 

40.5 

45.5 

77.9 

74.4 

62.1 

-115. 

-120. 

-83.0 

C28 

76 

025 

16 

09 

-148. 

5.45 

S 

72.51 

W 

49.0 


60.1 

87.2 

85. 


-101. 

-69.6 

C30 

76 

121 

11 

52 

21.7 

9.04 

S 

12.10 

E 

35.8 


36.4 

3.4 

-78.2 


165.9 

-71.3 

C31 

76 

137 

12 

47 

-8.2 

10.48 

S 

9.83 

E 

33.8 


37.1 

5.8 

-75.6 


169.9 

-80.1 

C35 

76 

319 

23 

16 

-216. 

19.37 

N 

89.01 

W 

68.2 


83.9 

106.6 

112.3 

• 

-71.4 

-46.8 

C36 

77 

256 

20 

34 

-13.7 

1.43 

S 

22.61 

1/ 

1.8 

5.6 



-153. 

-66.7 



HOI 

73 

171 

20 

24 

-160. 

20.01 

S 

72.41 

u 

50.7 

55.9 

87.0 

85.0 

66.5 

-112. 

-114. 

-84.0 

H02 

74 

192 


50 

-185. 

19.12 

N 

73.54 

E 

97.6 

92.2 

64.2 

63.7 

-109. 

71.0 

80.3 

88.9 

HO 3 

75 

3 

1 

46 

-237. 

28.47 

N 

92.45 

W 

73.3 

78.7 

82.7 

110.0 

121.3 

-59.7 

-61.5 

-37.1 

H04 

75 

12 

3 

17 

-161. 

59.59 

N 

53.65 

E 

86.1 

83.8 

47.8 

75.0 

-150. 

28.8 

31.6 

44.4 

H06 

73 

72 

8 

01 

-275. 

84.98 

S 

110.50 

W 

86.7 

86.6 

118.0 

84.0 

5.0 

-175. 

-175. 

-150. 

H07 

71 

192 

13 

28 

-201. 

20.16 

S 

72.16 

W 

50.5 

55.7 



66.2 

-112. 




128 





Table 18. (continued) 


Code Reference time Origin Epicenter (degree) Distances (degrees) Azimuths (degrees) 



yr:day 

hr:mn 

(sec) 

Latitude Longitude 

12 

14 

15 

16 

12 14 IS 16 

H08 

71 

140 

17 29 

-168, 

39.44 N 

18;62 W 

42.7 

43.1 



-175. -i.i 

H09 

71 

107 

07 05 

-247. 

49.45 N 

28.69 E 

69.0 

66.4 



-147. 30.8 

H10 

76 

66 

10 14 

-91.6 

49.22 N 

26.01 W 

52.3 

53.4 

32.6 

68.6 

177.8 -7.C -40.0 -2.3 

Hll 

76 

68 

14 43 

-46.0 

18.13 S 

11.83 W 

18.9 

15.5 

46.7 

28.1 

-36.2 159.5 -160. -86.4 

A01 

75 

113 

11 54 5) 


10.80 S 

31.30 W 

11.0 

15.5 

50,2 

46.1 

44.8-118, -133. -70.9 

A14 

73 

143 

18 58 


22.50 S 

35.30 W 




50.5 

-86.3 

A16 

72 

145 

09 24 


5.90 N 

6.80 E 




17.2 

-5.1 

A17 

73 

001 

11 33 


33.50 N 

33.10 W 

37.7 

40.0 

32.6 

62.7 

166,7 -20.5 -68.0-19,2 

A20 

74 

116 

09 20 


20.80 N 

27.00 W 

24.1 

26,2 

28.6 

51.3 

171.8 20.5 -94.1-28.5 

A28 

71 

288 

06 29 


7.00 N 

24.20 E 


42,9 



75.7 

A33 

72 

341 

13 30 


7.10 N 

1.05 E 




90.6 

108,4 

A34 

72 

138 

23 33 


7.50 N 

6.70 M 



21.0 

27.6 

-150. -28.5 

A42 

74 

166 

13 21 


23.10 N 

45.70 W 

34.0 

38.4 

44.7 

67.9 

141.4 -44.5 -82.7-35,0 

A44 

74 

084 

14 33 


61.70 N 

16.80 E 


70.4 


70.7 

16.5 26.1 

A08 

77 

107 

17 15 


26.30 N 

32.40 W 




48.5 

-91.8 


129 



130 


Table 18. (continued) 

^ Origin times and locations of C and H events were obtained by using 
arrival times listed In Table 19 and the velocity distributions given 
In Table 17. Source depths were assumed to be 9 km for X events *and 80 
km for H events. Locations of A events are -from Lamnleln et al. (1977) 
and were used for calculations of the distances and azimuths only. 

21 

' Alternate locations of those events which had arrival times observed at 
two stations only are as follows: 


Code 

Latitude (°) 

Longitude (°) 

C02 

7.91 S 

18.09 W 

C03 

36.30 S 

37.25 W 

H07 

24.06 N 

67.29 W 

H08 

45.38 S 

29.81 W 

H09 

62.61 S 

23.93 E 


3) 

' No distance and azimuth are given if there was no long-period seismogram 
available at a given station or if the seismic signal was not used in 
the seismic sections. 

^ Azimuth is the angle from the positive X axis of the seismometer to the 
direction toward the source from the station (see Table 1 for the 
orientation of the long-period seismometers). Positive is clockwise. 

^ Reference time of the seismic signal at station 14 Is 71:051; 11:54. 



Table 19. Observed and theoretical arrival times for events listed In Table 18 


1 ) 


Code Station 12 Station 14 Station 15 Station 16 



P obs 

P cal 

S obs 

S cal 

P obs 

P cal 

S obs 

^cal 

P obs 

P cal 

S obs 

S cal 

P obs 

■“cl S 

obs S cal 

CL2 

38.0 

37.6 

55.5 

52.8 












CL4 

51.5 

52.0 

67.5 

72.3 

44.5 

44.5 

60,0 

59.0 








CL5 









-0.5 

0.5 

24.5 

28.3 




CS3 

10.0 

10.8 

34.0 

34.8 












CS4 

31.0 

31.1 

56.0 

55.9 












CS5 

41.0 

40.6 

87.5 

85.0 

19.5 

20.0 


47.1 








CS6 

31.5 

31.5 

56.5 

56.5 

47.0 

47.0 

86.0 

81.0 

150.0 

147.0 

259.0 

258.8 




CS7 

39.0 

39.2 

79.0 

81.7 

14.3 

15.4 

44.0 

39.6 

134.5 

125.6 


232.8 

105.0 103.1 

193.4 

C01 

25.0 

25.5 

55,0 

57.6 

12.5 

12.5 

36.5 

36.6 

114.0 

114.0 

216.0 

214.3 

121.5 

122.4 

228.9 

C02 

19.5 

19.5 

50.0 

50.2 

6.5 

6.5 

28.0 

28.0 







• 

C03 

52.0 

51.7 

164.0 

164.0 

58.5 

58.5 

176.0 

176.0 








C04 

35.5 

35.5 

143.0 

142.9 

41.5 

42.4 


155.0 

15.0 

15.0 

107.0 

107.1 




C05 

114.5 

114.5 

• 

508.3 

123.5 

113.1 


505.1 

55.0 

55.0 


380.5 

89.0 

88.9 

451.0 

C06 


165.9 


293.2 

160.0 

160.0 


282,9 

51.0 

51.0 

88.0 

88.0 


185.6 

327.8 

C07 

90.0 

91.4 

286.0 

282.2 

94.5 

94.5 

292.0 

288.0 

7.0 

7.0 

132.0 

132.2 

131.5 

133.7 

363.2 


Table 19. (continued) 


Code Station 12 Station 14 Station 15 Station 16 



p 

K obs 

P cal 

S obs 

S cal 

P obs 

P cal 

S obs 

S cal 

P obs 

P cal 

S obs 

S cal 

P obs 

p «i 

S obs 

S c.l 

C16 

149.0 

155.1 

294.0 

293.1 

132.0 

133.2 

252.0 

254.8 

124.0 

124.0 

239.0 

238.7 

11.0 

11.0 

40.5 

40.5 

C23 

171.5 

175.8 

350.0 

348.8 

157.0 

155.2 

313.0 

312.4 

100.5 

100.5 

210.0 

216.8 

44.5 

44.5 

124.0 

118.9 

C25 

2.0 

4.1 


304.9 

14.5 

14.5 


325.8 

76.5 

76.5 


455.4 

53.5 

53.5 


405.8 

C26 

6.0 

1.5 

141.0 

126.6 

19.5 

20.0 

160.0 

158.9 

130.0 

130.0 

361.0 

361.0 

119.0 

119.0 


340.0 

C28 

50.0 

50.0 

196.0 

198.4 





154.0 

154.0 


90.7 

175.0 

175.0 


430.3 

C30 


171.7 

284.0 

284.0 






174.7 

286.0 

287.5 

46.5 

46.5 

64.5 

65.3 

C31 


134.4 

241.0 

241.0 






146.4 

262.0 

262.1 

27.5 

27.5 

54.0 

53.8 

C35 

49.0 

49.0 


251.4 

67.0 




98.0 

98.0 

345.0 

344.6 

158.5 

158.5 


467.8 

C36 

3.0 

3.0 

16.0 

15.8 

21.0 

21.0 


46.6 








267.1 

H01 


35.7 

182.0 

182.0 

54.0 

54.0 

215.0 

215.0 

156.0 

151.7 


399.7 


146.1 


388.7 

H02 


155.4 


434.4 

128.0 

141.2 


404.1 

57.0 

57.0 

242.0 

242.0 

55.5 

55.5 

240.0 

239.3 

H03 

34.0 

34.0 


245.7 

51.0 

50.6 

275.0 

276.5 

62.5 

62.4 


298.6 

133.0 

133.0 


442.3 

H04 


147.7 


393.4 

140.0 

140.9 


380.1 

24.0 

24.0 

161.0 

162.7 


114.7 

331.0 

331.0 

H06 

34.0 

35.3 


282.6 

35.0 

35.0 


282.1 

109.0 

111.3 


438.1 

27.5 

27.5 

266.0 

267.3 


H07 


-5.7 140.0 140.0 12.5 12.5 173.0 173.0 


U 

l\> 


fable 19. (continued) 


Code 


Station 12 


Station 14 


Station 15 


Station 16 


P obs P cal S obs S cal P obs P ca'i ^obs ^cal P obs P cal ^obs S cal P obs P cal S obs S cal 


H08 

-1.1 123.5 123.5 

0.0 

o 

• 

o 

126.0 

126.0 




H09 

10.0 10.0 210.0 210.0 

1.5 

1.5 

192.0 

192.0 




H10 

109.7 262.5 260.6 

113.5 

113.5 

267.5 

267.4 

38.5 

38.5 132.5 135.9 

164.5 

Hll 

32.8 91.7 

20.0 

20.0 

69.0 

69.4 

139.5 

135.3 271.0 271.0 

66.0 67.1 150.0 

A01 

106. 



116. 



213. 

192. 

A14 








22.5 

A16 








170. 

A17 

29. 



29. 



50. 

107.5 

A20 

27. 



30. 



37. 

105.5 

A28 




38. 





A33 








218.5 

A34 







39. 

52.5 

A42 

58.5 



71. 



93. 

170.5 

A44 




26.5 




51. 

A08 








4. 


363.3 


U 



134 


Table 19. (continued) 

^ Reference times are listed In Table 17. Arrival times are given 
for those signals only which were used In the seismic sections. 

The seismic records of C and H events are aligned In the sections 
according to either the observed or, If underlined, the calculated 
arrival times of the congressional waves. The seismic records of 
A events are aligned according to the observed shear-wave arrivals. 
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Table 20. Arrival times observed by other Investigators.^ 


a) Artificial Impacts 


Event Sta. Latham etaU Toksozetal. Koyama 4 Naka- Jarosch (1977) 


Code 

(1972b) 

(1974) 

mura (1979) 2) 





P 

S 

P S 

P 

S 

P 

S 

CL4 

12 

-0.8 

3.2 

-2.0 


3.4 




14 

-1.0 

-3.3 

-0.3 

0.3 


0.5 

-2.4 

CL5 

15 

-0.5 


-0.5 


-9.7 


-10.3 

CS3 

12 

0.4 

-2.5 

-1.0 

0.0 

-3.3 

-0.2 


CS4 

12 

0.1 

-4.9 


0.2 

-3.1 

0.2 

1.0 

CS5 

12 

-2.6 


-2.6 

1.0 

-2.1 

4.6 

1.4 


14 

0.0 

43.4* 

-0.5 

0.1 

44.9* 

1.1 

* 

37.5 

CS7 

12 




0.2 

2.2 

3.4 

-3.3 


14 




0.5 

•6.2 

0.1 

-16.7 


15 







210.9* 


16 






-0.8 

185.2* 

b) 

Meterold Impacts 


• 




Event Sta. 

Nakamura et al. 

Goins (1978) 

Koyama & Naka- 





(1977) 


mura 

(1979) 2 * 





P 

S 

P S 

p 

S 



C01 

12 

0.2 


0.2 7.7 

0.3 





14 

-0.2 

0.5 

0.0 0.3 

-0.1 

0.5 




15 

0.1 

4.0 

0.3 1.0 

0.4 

1.0 




16 -1.3 220* 


0.9 


-1.1 222.5 





136 


Table 20. (continued) 


C 05 

12 

- 0.5 


0.5 


- 0.5 



14 

0.4 

630 * 

0.3 


0.5 

630 * 


15 

4.7 

380 * 

- 8.7 


- 8.5 

380 * 


16 

- 0.2 

530 * 

• 12.3 


- 12.9 

530 * 

C 06 

12 

19 c . 1 * 


196 . 4 * 


196 . 1 * 

338 * 


14 

16.4 

280 * 

18.1 


17.5 

307 * 


15 

. 1.9 

2.0 

0.3 

7.5 

2.0 

2.0 


16 

192 * 

320 * 

199 . 5 * 


193 * 

* 

331 

C 07 

12 



- 2.4 





14 

0.1 

8.0 

- 0.2 

- 1.8 




15 

15.1 

- 2.0 

14.3 





16 

- 1.1 

360 * 

- 0.2 




C 16 

12 

13.1 







14 

9.8 

18.0 






15 

0.0 

1.0 






16 

0.1 

9.5 





C 23 

12 

0.1 


0.3 

- 58.0 

0.2 

2.0 


14 


- 13.0 

1.0 


1.0 

7.0 


15 

0.3 

10.0 

0.3 


0.3 

10.0 


16 

- 0.2 

- 4.0 

0.1 

6.5 

- 0.1 

- 4.0 

C 25 

12 

0.3 


- 0.7 


0.3 

282 * 


14 

0.4 

320 * 

1.0 


1.6 

320 * 


15 

1.2 

510 * 

1.0 

^ 

1.3 

* 

416 


16 

0.3 

460 * 

0.1 


0.3 

375 * 
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Table 20. (continued) 


C26 

12 



0.9 

0.6 

37.0 


* 


14 



-1.2 

-0.5 

38.5 




15 



-2.1 

-1.8 





"16 



7.5 

0.5 




C28 

12 



1.1 -2.5 






15 



0.5 372.2* 






16 



-4.3 





c) 

Shal iow moonquakes 






Event Sta. 

Nakamura et al. 

Goins (1978) 

Koyama & Naka- 

Lamm! el n 



(1979) 2) 


mura (1979) 2) 

(1977) 3) 



P 

S 

p S 

P 

s 

P 

S 

HOI 

12 

56.0* 

-2.0 

55.0* 3.3 

56.0* 

5.5 

0.0 



14 

13.0 

4.0 

12.5 

13.0 

-7.5 

17.0 

165.3 


15 

148* 

360* 

145.7* 

147.5* 

346* 

98.0* 

316.5 


16 

191* 

430* 

198.5* 412.5* 


422.5* 



H02 

12 

139* 

459* 

138.5* 

138.7* 

* 

417. 

167.2 



14 

1.0 

439* 

-2.7 

0.6 

* 

407 

153.4 



15 

1.0 

13.0 

-0.5 

0.5 

21.0 

85.0 

184.5 


16 

4.5 

30.0 

-5.2 

4.5 

30.0 

81.8 

179.5 

H03 

12 

0.0 


-0.4 269* 

-0.5 

256* 

0.0 

217.2 


14 

0.0 

5.0 

0.3 

0.2 

10.0 

17.8 



15 

-1.5 

300* 

-2.0 

-1.4 

★ 

315 

29.0 



16 

2.0 

* 

510 

-5.5 453* 

2.0 

510* 

100.6 
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Table 20. (continued) 


H04 

14 

-2.0 

368* 



•2.5 

356* 




15 

-4.0 

1.0 



6.0 

1.0 

0.0 

130.9 


16 

125* 

0.0 



125* 

0.0 

94.4 

295.2 

H06 

12 

0.0 

267* 

0,1 

* 

272 

0.2 

267* 

6.5 

239.3 


14 

-1.0 

280* 

0.9 


0.5 

272.5* 




15 

-9.0 

465* 

-9.3 


-9.0 

398.5* 

71.0 



16 

-0.5 

4.0 

0.3 

-6.6 

0.5 

2.6 

0.0 

230.9 

H07 

12 

-2.0* 

0.0 






126.8 


14 

10.0 

-8.0 





0.0 

151.3 

H08 

12 

0.0* 

-0.5 






123.3 


14 

0.0 

-2.0 





0.0 

123.8 

H09 

12 

0.0 

2.0 





7.0 

207.7 


14 

0.5 

-10.0 





0.0 

192.2 

H10 

12 

* 

118 

0.5 

110,8* 

-0.5 

110.8* 

0.5 




14 

-0.5 

0.5 

-0.2 

1.2 

-0.5 

0.5 




15 

-0.5 

-8.5 

0.7 

3.3 

-0.5 

-2.5 




16 

168* 

347* 


346* 

* 

161 

355* 



Hll 

12 



98.8* 







14 

2.0 

0.0 

4.2 

1.3 

2.0 

0.0 




15 

10.5 

-5.0 

10.0 

-5.9 

9.5 

-5.0 




16 

-6.0 

-16,0 

3.9 


-6.0 

-16.0 




^Differences between arrival tines read by other investigators and used 
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Table 20. (continued) 

in this study (see Table 19) are given in seconds. Arrival times 
read only by other investigators are marked with * and were 
measured from the reference times given in Table 18. 

21 

Observed arrival times were used in the publications but were un- 
published; the arrival times were provided by the first authors. 
31 

'Only relative arrival times were given hy Lammlein (1977). 
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Figure 1. Compressed time-scale seismic records for the three 
major types of natural seismic sources (from Nakamura et a!. , 
1974 ). 


preceding page blank not 
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Figure 2. Simplified block diagram of the feedback controlled, 
direct digitizing long-period lunar seismographs. 


FIGURE 2. 


Calibration 



Tidal output 



145 


Figure 3. Average of the positive and the sign-reversed 
negative calibration pulses in the peaked mode. The 
maxima of the calibration pulses are given in Table 3. 
The minima are 0. 




FIGURE 3 
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Figure 4. Average of the positive and the sign-reversed 
negative calibration pulses In the flat mode. See Figure 3 
for explanation; the minima are also given In Table 3. 
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{ 


FIGURE 4. 
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Figure 5. Schematic flow chart for calculating the seismograph 
parameters. 
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FI gum 6. Real and Imaginary parts of the partial derivatives 
of the transfer function to a step of acceleration In the 
peaked mode. The parameters, a^, and their units are 
listed In Table 2. The calculations were carried out with 
seismograph constants obtained for component 15X (s>e 
Table 4) and with KK^K^IO DU/cm. The minima and maxima 

3 

are given In DUsec /cm/unit of each parameter. 
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Figure 7. Real and Imaginary parts of the partial derivatives 
of the transfer function to a step of acceleration In the 
flat mode. See Figure 6 for explanation. 
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Figure 8. Partial derivatives of the response to a step of 

acceleration in the peaked mode. See Figure 6 for explana- 

2 

tion. The minima and maxima are given in DUsec /cm/unit 
of each parameter. 
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Figure 9. Partial derivatives of the response to a step of 
acceleration in the flat mode. See Figure 8 for 
explanation. 
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Figure 10. Partial derivatives of the amplitude and phase 
response to a step of acceleration in the peaked mode. 

See Figure 6 for explanation. The minima and maxima are 

3 

given in DUsec cm/unit of each parameter for the amplitude 
spectra and in radian/unit of each parameter for the phase 
spectra. 
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Figure 11. Partial derivatives of the amplitude and phase 
responses to a step of acceleration In the flat mode. 
See Figure 10 for explanation. 



Figure 12. Observed and theoretical amplitude and phase 
responses to a step of acceleration for components 15Z 
15X In the peaked mode. The minima and maxima are - 
and , respectively, for the phase responses. 



Figure 13. Absolute values of the differences between the 
observed and theoretical amplitude responses to a step 
of acceleration for components 15X and 15Z, and the rela- 
tive errors In relation to the amplitude responses In the 
peaked mode. The relative errors are clipped In order to 
utilize the full amplitude scale In the frequency band 
of Interest. 
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* • t 
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FREQUENCY (HZ) 



I5X 

FIGURE 13. 


I5Z 



Figure 14. Amplitude responses to an impulse of displacement 
in the peaked mode. The calculations were done with con- 
stants listed in Table 4 and with KK^K^’l DU/cm. The units 
are DUsec/cm. 








Figure 15. Absolute values of the relative differences 
between the amplitude responses of Figure 12 and the 
amplitude response of component 15Z in the peaked mode. 
The differences were not calculated at frequencies where 
the amplitude response is less than 3 % of the peak 
amplitude response. The amplitude responses of the 
various components were normalized by multiplying them 
with scale factors given in Table 6. 
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Figure 16. Absolute values of the relative differences 
between the amplitude responses of Figure 12 and the 
amplitude responses of component 15X in the peaked mode. 
See Figure 15 for explanation. The scale factors are 
given in Table 7. 


DEVIATION IN AMPLITUDE RESPONSE RELATIVE TO I5X 

(PEAKED MODE) 




FREQUENCY (Hz) 


FIGURE 16. 



Figure 17. Schematic flow chart for calculating the 
waterlevel parameters, the reliable frequency bands, 
and the relative errors between the observed and 
theoretically amplitude responses to a step of 


acceleration. 




FIGURE 17 





















Figure 18. Schematic flow chart for calculating the transfer 
functions of the seismographs and for adjusting the 
transfer functions with a given waterlevel parameter. 

Figure 19. Schematic flow chart for performing wavelet 
deconvolution by spectral division. 



176 
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FIGURE 18. 


FIGURE 19. 




























Figure 20. Schematic flow chart for stacking the amplitude 
spectra of the seismic records corrected for Instrumental 
response in the peaked mode. 

Figure 21. Schematic flow chart for calculating the total 
transfer function of the instrumental and near-surface 


effects. 
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FIGURE 2. . 
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Figure 22. Illustration of the noise of the analcg-to-digital 
converter at station 14. 
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Figure 23. Sum of the amplitude spectra for the horizontal 
components (X and Y) of station 12. The selected events are 
listed in Tables 8 and 9. The following abbreviations are 
used: A - deep moonquakes, H - shallow moonquakes, 

C - impacts (artificial and meteoroid impacts). The minima 
and maxima of the curves are given in narantheses. Units are 


DUsec. 



12Y, C C. 3/53, 16.27) 
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FIGURE 23. 
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Figure 24. Sum of the amplitude spectra 
components of station 14, See Figure 


for the horizontal 
23 for explanation. 



UY, C (.2866,35.06) 
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Figure 25. Sum of the amplitude spectra for the horizontal 
components of station 15., See Figure 23 for explanation. 
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Figure 26. Sum of the amplitude spectra for the horizontal 
components of station 16. See Figure 23 for explanation. 



16Y, C (,08112,11.32) 
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FIGURE 26. 
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Figure 27. Sum of the amplitude spectra for the vertical 
components (Z) of stations 12 and 14. See Figure 23 for 
expl anation. 
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Figure 28. Sun of the anplitude spectra ror the vertical 
conponents of stations 15 and 16. See Figure 23 for 


explanation. 


16Z , C 
. 1969,2.99) 
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Figure 29. Sun of the anplitude spectra for station 12. HOR 

is the horizontal component calculated with equation (28). 

The amplitudes are smoothed with moving averages in a 
rectangular winaow of 0,055 Hz. 



FIGURfc 29. AMPLITUDE SPECTRA AT STATION 12 
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Figure 30. Sum of the amplitude spectra for station 14, See 
Fiaure 29 for explanation. 



FIGURE 30. AMPLITUDE SPECTRA AT STATION 14 
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Figure 31. Sum of the amplitude spectra for station 15. See 
Figure 29 for explanation. 
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Figure 32. Sum of the amplitude spectra for station 16. See 
Figure 29 for explanation. 
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Figure 33. Two examples for instrumental and whitening deconvo- 
lution. The P wave arrival times at which the records are 
aligned are given in Table 19. 



FIGURE 33 
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Figure 34. 
(X to Y 


Ratio of the average horizontal spectral amplitudes 
. See Figure 29 for explanation. 
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Figure 35. Ratios of the horizontal and vertical spectral ampli- 
tudes at station 12. The ratios were calculated individually 
and were averaged for deep moonquakes (A), shallow moonquakes 
(H), and impacts (C). Curve 'Total' is the spectral ratio 
curve calculated from the average spectral amplitudes (see 
Figures 29 through 32). The amplitudes are clipped at long 
periods in order to utilize the full amplitude scale in the 
frequency band of interest. The average ratios were smoothed 
with moving averages in a window of 0.55 Hz. The dot-dash 
and the dashed curves are the theoretical curves calculated 
for models 'A' and '9', given in Tables 12 and 13, 


respectively. 
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Figure 36. Ratios of the horizontal and vertical spectral ampli- 
tudes at station 14, See Figure 34 for explanation. 
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Figure 37. Ratios of the horizontal and vertical spectral ampli- 
tudes at station 15. See Figure 34 for explanation., 
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Figure 38. Ratios of the horizontal and vertical spectral ampli- 
tudes at station 16. See Figure 34 for explanation. 
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Figure 39. Near-surface shear velocity profiles obtained 
for the Apollo 12, 14, 15, and 16 landing sites. The 
model parameters are given in Table 13. The center of 
the layers is marked with a cross. The dashed line is 
the piecewise linear shear velocity function described 
in Table 14. 
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Figure 40. Travel time curves of the refracted congressional 
waves for station 12. The model parameters are given in 
Table 13. The first arrival from the Lunar Module lift-off 
is indicated with a star. 

Figure 41. Travel time curves of the refracted congressional 
waves for station 14. See Figure 39 for explanation. 



Travel time (seconds) Travel time (seconds) 



FIGURE 41. 
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Figure 42. Travel time curves of the refracted compressional 
waves for station 15. See Figure 39 for explanation, 

Figure 43. Travel time curves of the refracted compressional 
waves for station 16. The model parameters are given in 
Table 13. 



T ravel tine (seconds) Travel tine (seconds) 


FIGURE 42. 




FIGURE 43 
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Figure 44. Velocity nodels used in this study for the whole 
noon. The velocities of the model with thicker crust and 
higher mantle velocities (continuous line) are given in 
Table 16. The velocities of the model with thinner crust 
and lower mantle velocities (dashed line) are given in 
Table 17. 






Figure 45. Velocity models used in this study fc the lunar 


crust and upper mantle only. See Figure 43 for evp' mation. 




FIGURE 45. 


DEP 
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Figure 46. Locations of impacts as determined in this study. The 
base nap is the whole noon in an equal area projection. Table 
17 lists these events. The location of the seismic stations 
are also given. 


2 24 



FICURE 4b. 




225 


Figure 47. Epicenters of shallow moonouakes as determined in 
this study. See Figure 45 for explanation* 
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FIGURE 47. 



SEISMIC SECTION CAPTIONS 


IMPACTS (R, 0.4-1 .5HZ) ' Radial component of ground motion, distance 
range is 0 to 50". Seismograms were corrected for instrument response 
and near-surface effects, filtered with a 3-pole Butterworth filter, 
and rotated into the above direction with azimuths given in Table 18. 
The distances, at which the records are aliened, are given in Table 18. 
and are indicated by drawing a line from the given distance to the be- 
ginning of the seismic record. The arrival times of the direct waves, 
congressional (P) or shear (S), at which the seismic records are lined 
up, are given in Table 19. The seismic records are identified with 
a sequence numoer, with a four-character word (first tnree characters 
are the codes given in Tables 9 and 13, the forth character is the 
second digit of the station number), and with a number giving the 
scale factor which was used to normalize the amDlitudes by division. 
Glitches were zeroed. If a component was not operating normally (see 
Table 9), then the seismogram was zeroed. The part of the seismic 
records where the amplitudes are clipped were not plotted. The 
following arrival times were calculated with velocity distributions 
given in Table 17 and in Figures 43 and 44 and were drawn with contin- 
uous lines: S - direct shear wave arrival; PS-P to S reflection at 

the free surface; SS, 3S and 4S - shear wave is reflected once, 
twice and three times at the free surface; S25-peg leq multiple in the 

upper crust; S45-oeg leg multiple in the whole crust; SLC - oeg leg 

?00 

multiple in tne lower crust; $4: P - S to P conversion at 200 km. 

IMPACTS ( T,0.4-l . 5HZ) 1 Transverse component of ground motion, distance 
range is C to 50°. See Seismic Section 'IMPACTS (R, 0.4-1 .SHI ) 1 
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for explanation. 

'IMPACTS (Z, 0.4-1 . 5HZ ) ' Vertical component of ground motion, distance 
range is 0 to 50°. See Seismic Section 'IMPACTS ( R ,0.4-1 . 5HZ) ' 
for explanation. 

'IMPACTS (R, 0.4-1 . 5HZ ) ' Radial component of ground motion, distance 
range is 60 to 160°. See Seismic Section 'IMPACTS (R,0.4-l . 5HZ) ' 
for explanation. 

'IMPACTS (T, 0.4-1 .5HZ) ' Transverse component of ground motion, distance 
range is 60 to 160°. See Seismic Record Section 'IMPACTS (R, 0.4-1 .5HZ) ' 
for explanation. 

'IMPACTS (Z, 0.4-1 . 5HZ) ' Vertical component of ground motion, distance 
range is 60 to 160°. See Seismic Section 'IMPACTS (R, 0.4-1 . 5HZ) ' 
for explanation. 

'SHALLOW M00NQUAKES (R, 0.4-1 . 5HZ ) ' Radial component of ground motion; 

distance range is 0 to 120°. See Seismic Section 'IMPACTS (R, 0.4-1 .5HZ) ' 
for explanation. 

'SHALLOW M00NQUAKES (Z, 0.4-1 . 5HZ ) ' Vertical component of ground motion, 
distance range is 0 to 120°. See Seismic Section 'IMPACTS ( R ,0 . 4-1 . 5HZ) ' 
for explanation. 

'DEEP M00NQUAKES (R, 0.4-1 . 5HZ ) ' Radial component of ground motion, 

distance range is 0 to 100^*. See Seismic Section 'IMPACTS (R, 0.4-1 . 5HZ) ' 
for explanation. 

'DEEP M00NQUAKES (T, 0.4-1 . 5HZ) ' Transverse component of ground motion, 
distance range is 0 to 100°. See Seismic Section 'IMPACTS (R, 0.4-1 . 5HZ ) ' 
for explanation. 
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'DEEP MOONOUAKES ( Z) ' Vertical component of ground motion as 

recorded with seismographs operated in the peaked mode, distance range 
is 0 to 100°. See Seismic Section 'IMPACTS (R, 0.4-1 .5HZ) ' for 
explanation. 
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